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Chapter  I 
INTRODUCTION 


1.1  Turbulence-Induced  Pressure  Helds 

The  pressure  induced  by  boundary- layer  turbulence  has  generated  much 
concern,  arising  primarily  from  the  unsteady  loads  it  induces  on  flexible 
panels  and  the  noise  that  it  radiates  (due  to  pressure  fluctuations  in  the 
fluid,  or  indirectly  from  panel  vibration).  Panel  vibration  is  driven  by  the 
spatially  integrated  turbulent  boundary  layer  (TBL)  pressure  field,  weighted 
against  the  structure's  response  function,  so  that  the  integral  scales  of  the 
pressure  fluctuations  are  important  measures  of  boundary  layer  turbulence. 
Essentially,  the  correlation  area  of  the  pressure  producing  turbulent  eddies 
determines  the  mean-square  level  of  the  applied  force  field. 

The  subject  of  wall  pressure  fluctuations  beneath  turbulent  boundary 
layers  has  been  extensively  studied  (for  reviews  see  Corcos  (1964),  Haddle 
and  Skudrzyk  (1969),  and  Willmarth  (1975)1.  From  these  analytical  and 
experimental  investigations,  it  is  now  widely  accepted  that  the  energy 
associated  with  TBL  wall  pressure  fluctuations  is  distributed  over  wide 
frequency  and  wavenumber  ranges.  Three  distinct  wavenumber  regions  of 
the  turbulent  pressure  fluctuations  have  been  identified;  they  are,  in  order 
of  increasing  wavenumber,  the  sonic,  low  wavenumber  and  convective 
regions,  respectfully. 

The  lowest  wavenumber  region  (longest  wavelengths)  of  the  TBL 
pressure  fluctuations  consists  of  wavenumbers  less  than  and  equal  to  the 
sonic  wavenumber  (k-u/c)  of  the  fluid,  where  a)  is  the  radian  frequency  and 
c,  the  sound  speed  in  the  fluid.  This  pressure  field,  although  very  weak, 
radiates  as  true  sound  and  is  caused  by  the  interaction  of  the  turbulence 
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with  itself.  Many  investigators  have  addressed  this  sound  field  with  major 
contributions  from  Lighthill  (1952,  1954),  Haddle  and  Skudrzyk  (1969).  and 
Tam  (1975).  A  few  of  the  many  other  experimental  studies  of  this 
TBl-induced  sound  field  include  Nishi  et  at.  (1970),  Skudryzk  and  Haddle 
(1960),  Meecham  (1965),  Vecchio  and  Wiley  (1973)  and  Lauchle  (1976). 

Wall  pressure  fluctuations  at  wavenumbers  greater  than  sonic 
wavenumber  (k)  and  smaller  than  the  convective  wavenumber  (kc=o>/uc, 
where  Uc  is  the  convection  velocity),  describe  the  low  wavenumber  region 
of  the  TBL  pressure  fluctuations.  It  is  presumed  that  large-scale  coherent 
motions  in  the  TBL  are  characterized  by  bursts  and  sweeps  of  turbulent 
energy  and  contribute  significantly  to  the  pressure  fluctuations  in  this  low 
wavenumber  region  (eg..  Kline  etal.  (1963),  and  Blackweider  and  Kaplan 
(1976)).  Although  comprising  a  very  small  portion  of  the  total  TBl's 
wavenumber -frequency  spectral  energy  it  is  thought  that  these  large-scale 
motions  in  the  TBL  are  major  causes  of  excessive  structural  fatigue 
Analytical  models  for  the  pressure  field  at  these  low  wavenumbers  are 
given  by  Chase  (1980),  Maestrello  (1965,1967),  and  Wetting  (1976). 
Extensive  measurements  of  the  wall  pressure  fluctuations  in  this  low 
wavenumber  region  do  not  exist  because  of  the  inherent  experimental 
difficulty  with  such  measurements.  For  example,  spatial  filters  are  needed 
to  eliminate  both  very  low  and  high  convective  wave  contributions. 
Farabee  and  Geib  (1976),  Willis  (1970),  Blake  and  Chase  (1971),  and  Martin 
and  Leehey  (1977)  have  performed  some  of  these  measurements;  however, 
reasonable  agreement  between  experimental  results  and  theoretical 
predictions  are  less  than  satisfactory. 

The  final  region  of  the  TBL  wavenumber  spectrum  consists  of 
wavenumbers  that  lie  in  a  narrow  ridge,  centered  about  the  convective 


wavenumber  (kc)  of  the  flow.  Pressure  fluctuations  in  this  convective 
region  arise  from  the  interaction  of  the  turbulence  and  the  mean  shear. 
Current  models  suggest  that  this  region  is  composed  of  the  dominant  TBL 
fluctuating  energy  and  comprises  the  major  part  of  the  r.m.s.  wall 
pressure.  This  so-called  'convective  ridge'  energy  arises  because 
small-scale  turbulent  motions  within  the  TBL  convect  over  a  surface  at  a 
preferred  convection  velocity  (typically  0.7  to  0.8  times  the  free-stream 
velocity).  However,  these  small-scale  eddies  are  correlated  for  very  short 
distances  and  decay  due  to  the  viscous  damping  effects  of  the  fluid,  in 
fact,  the  distances  that  these  pressure  producing  eddies  convect  is  directly 
proportional  to  their  size.  The  pressure  field  along  this  convective  ridge 
has  perhaps  received  the  most  attentioa  beginning  with  the  theoretical 
modeling  by  Kraichnan  (1956a.  b).  Wlllmarth  (1975)  gives  an  account  of  tne 
most  significant  contributions  to  this  area. 

Although  the  energy  contained  in  this  pressure  ridge  dominates  the 
wavenumber-frequency  spectrum,  its  accurate  measurement  over  a  wide 
frequency  range  is  not  a  trivial  matter.  Early  experimental  studies, 
particularly  those  of  Willmarth  (1956)  and  Corcos  (1963,1964)  recognized 
the  necessity  of  'small'  pressure  transducers  in  conducting  TBL  wall 
pressure  fluctuation  measurements.  The  convective  domain  is  primarily 
dominated  (spatially)  by  small-scale  eddies.  Therefore,  finite-sizeu 
pressure  sensors— which  average  over  space-need  to  be  at  least  as  small 
as  the  smallest  eddies;  otherwise,  an  attenuation  of  the  high-frequency 
(small  wavelength)  spectral  contributions  (particularly  those  in  this 
convective  ridge)  will  occur.  Kolmogoroff  (see  Hinze,  1959,  for 
discussion)  deduced  both  length  and  velocity  scales  for  these  smallest 
eddies.  He  found  that  the  size  of  the  smallest  eddies  is  inversely 


proportional  to  the  l/4th  power  of  the  dissipation  rate  of  turbulent 
energy.  Conversely,  he  found  that  the  velocity  scale  is  directly 
proportional  to  the  l/4th  power  of  this  same  dissipation  rate.  Therefore, 
by  combining  these  two  quantities  in  a  Reynolds  number  fashion  it  is  easily 
demonstrated  that  these  smallest  eddies  scale  as  u/ur,  where  o  is  the 
kinematic  viscosity  and  uT  is  the  friction  (or  shear)  velocity.  Emmerling 
(1973)  was  the  first  investigator  to  verify  experimentally  that  the 
smallest  scales  of  the  turbulence  are  indeed  the  same  order  as  the  viscous 
length  scale  0/1%.  The  parameter  that  most  accurately  describes  the 
degree  of  spectral  resolution  for  a  given  wall  pressure  measurement  is 
therefore,  d+=urd/o,  where  d  is  the  pressure  transducer's  active  diameter. 
Clearly,  d*  must  be  of  order  one  to  resolve  the  smallest  scale  pressure 
fluctuations,  e.g.  Emmerling  (1973). 

As  stated,  many  of  the  early  wall  pressure  fluctuation  measurements 
(Willmarth  (1956),  Willmarth  and  Wooldridge  (1962),  and  Bull  and  Willis 
(I960,  among  others)  suffered  from  transducer  spatial  averaging  problems. 
In  order  to  alleviate  these  problems,  Corcos  (1963)  developed  a  correction 
factor  that  could  be  applied  directly  to  the  measured  TBL  spectral  levels. 
Corcos'  correction  factors,  however,  were  based  on  the  similarity 
parameter  <or/uc  (where  r=d/2)  and  not  on  d*.  In  another  investigation 
Blake  (1970)  customized  an  l/8-inch  condenser  microphone  by  placing  a 
cap,  perforated  with  an  1/32-inch  diameter  hole  over  its  center.  With  this 
type  of  approach,  Blake  hoped  to  resolve  the  high-frequency  part  of  the  TBL 
wall  pressure  spectrum.  Blake’s  results  clearly  indicated  that  the  TBL 
pressure  fluctuations  contained  a  larger  frequency  content  (smaller  scale 
pressure  fluctuations)  than  what  had  been  measured  previously.  However, 
the  effect  of  the  hole  on  the  TBL  wall  pressure  spectrum  was  not  clearly 
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understood.  That  is,  was  the  observed  higher  frequency  content  of  the  TBL 
wall  pressure  spectrum  attributed  to  the  pressure  fluctuations  in  the  TBL 
or,  was  it  due  to  some  other  anomalous  contribution  from  the  pin-hole. 
Bull  and  Thomas  (1976)  later  addressed  this  question  by  comparing  the  TBL 
wall  pressure  spectrum  measured  by  a  microphone  with  a  pin-hole  cap 
(pin-hole  microphone),  to  that  measured  with  a  quartz-type  transducer  of 
the  same  diameter  as  the  pin-hole.  Their  findings  indicated  that  the 
increase  in  the  frequency  content  of  the  TBL  wall  pressure  spectrum  were, 
in  part,  attributable  to  the  hole  and  the  small  resonance  cavity  that  existed 
between  the  perforated  cap  and  transducer  diaphragm.  Subsequent 
investigations  of  the  TBL  wall  pressure  spectrum  have  followed,  together 
with  controversies  centered  on  the  parameters  tnat  best  describe  the 
ultimate  in  transducer  resolution  and  the  scaling  of  the  TBL  wall  pressure 
spectrum. 

Without  the  use  of  pin-hole  type  pressure  transducers,  Schewe  (1903) 
was  the  first  investigator  to  obtain  wall  pressure  fluctuation  data  for 
d*<50.  His  lowest  quoted  value  is  a  d*  of  19,  and  he  states  that  his 
measurements  of  the  TBL  wall  pressure  spectrum  are  representative  of  a 
true  point  measurement;  that  is,  there  should  be  no  attenuation  of  the  TBL 
spectrum's  high-frequency  content.  In  this  investigation,  TBL  wall 
pressure  spectra  are  obtained  using  a  unique  facility  and  data  processing 
technique.  The  transducer  spatial  resolution  achieved  is  significantly 
better  than  any  previous  investigation;  a  d+  of  2.1  is  demonstrated. 

1.2  Statement  of,  the.  Objective 

The  uniqueness  of  the  TBL  wall  pressure  measurements  conducted  in  this 
investigation  lie  in  the  fluid  and  the  facility  in  which  they  were  made.  All 
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data  obtained  in  previous  investigations  were  obtained  in  either  water  or 
air.  This  investigation,  however,  uses  glycerine  as  the  working  fluid. 
Glycerine  has  a  kinematic  viscosity  (u),  at  room  temperature,  of  almost 
1000  times  that  of  water  and  roughly  50  times  that  of  air.  Because  of  this, 
the  smallest  normalized  transducer  diameter  (di'=2.l)  to  date  has  been 
achieved.  Therefore,  the  high-freguency  spectral  energy  of  the  TBL  can  be 
investigated  by  the  acquisition  and  comparison  of  the  spectral  data  at  a  d* 
which  is  a  factor  of  nine  (9)  smaller  than  that  achieved  in  the  Schewe 
investigatioa  Also,  the  unique  way  in  which  the  flow  Reynolds  number  is 
altered  (viscosity  is  primarily  the  independent  variable  rather  than  the 
mean  flow  velocity)  allows  for  a  better  understanding  of  the  scaling  laws 
that  characterize  the  TBL  wall  pressure  spectrum. 

Investigations  involving  single  point  measurements  of  the  TBL  wall 
pressure  spectra  are  often  plagued  with  low-frequency  facility  noise.  The 
measurement  technique  used  in  this  investigation  minimizes  these  effects 
while  maintaining  the  integrity  of  the  measured  TBL  pressure  spectrum.  In 
addition,  a  measurement  involving  the  coherence  function  between  a  set  of 
(co-planar)  transducer  difference  signals  is  shown  to  validate  the  TBL  wall 
pressure  spectra  (especially  at  low  frequencies). 

1.3  Scope  of  the  Thesis 

This  chapter  introduced  some  fundamental  aspects  of  turbulence  and  the 
basic  types  of  noise  associated  with  it.  The  problems  to  be  addressed  in 
subsequent  chapters  were  also  described  and  the  anticipated  results 
discussed. 

Chapter  2  concentrates  on  the  techniques  used  to  measure  the  turbulent 
boundary  layer  (TBL)  wall  pressure  fluctuations.  This  selected  technique  is 


an  improvement  on  a  method  previously  used  by  Horne  and  Hansen  (1981),  in 
which  an  axisymmetric,  flush-mounted,  pressure  transducer  array  is  used 
to  minimize  the  farfield  acoustic  effects  in  the  TBl  wall  pressure 
spectrum. 

The  types  of  pressure  transducers  used  through  this  investigation  are 
described  in  Chapter  3.  Also  described  is  transducer  calibration  scheme. 
The  measured  sensitivities  and  relative  phase  responses  of  the  pressure 
transducers  are  given. 

Chapter  4  presents  data  on  the  measured  TBL  wall  pressure  spectrum  at 
a  number  of  flow  Reynolds  numbers  and  compares  these  results,  in 
non-dimensional  form,  to  those  obtained  by  other  investigators.  Also 
described  are  important  scaling  parameters  that  aid  in  the  understanding  of 
the  TBL  wall  pressure  spectrum.  This  chapter  concludes  with  a  discussion 
of  the  facility  generated  background  noise. 


Chapter  2 

MEASUREMENT  TECHNIQUE  AND  ANALYSIS 
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Measurements  of  turbulent  boundary  layer  (TBL)  pressure  fluctuations 
with  flush  (wall)  mounted  pressure  sensors  are  frequently  contaminated 
with  facility  generated  and  vibrationally  induced  noise.  A  novel  technique 
described  by  Wambsganss  (1979)  utilized  a  two  sensor  arrangement  that 
discriminated  between  the  propagating  (facility)  and  nonpropagating  (TBL) 
pressure  components  in  a  fully,  turbulent  pipe  flow.  However,  this  early 
technique  involved  uncertainties  when  combining  various  source  (TBL) 
terms. 

Refined  techniques  were  developed  later  by  Wilson  (1979),  and 
investigators  Horne  and  Hansen  (1981).  Both  of  these  techniques  involved 
an  axisymmetric,  three  sensor  configuration  that  resolved  earlier 
ambiguities  and  greatly  improved  the  integrity  (Wilson)  and  moreover,  the 
ease  of  acquiring  data  (Horne  and  Hansen).  All  of  these  methods  however, 
neglected  the  vibrationally  induced  pressure  signal  that  can  dominate  the 
output  of  flush-mounted  pressure  sensors  (especially  at  low  frequencies). 
The  following  describes  a  more  complete  analysis  and  represents  a  further 
refinement  to  the  latest  method  of  Horne  and  Hansen.  This  new  refinement 
includes  the  effect  of  tunnel  wall  vibration  and  offers  a  correction,  if 
needed-due  to  vibrational  effects— to  the  measured  TBL  pressure  spectra. 


Consider  the  flush-mounted  array  depicted  below  in  Figure  2.1.  The 
outputs  a(rj ,  e,  t),  b(rj ,  e,  t),  c(r| ,  e  ,t)  define  the  time  varying  output 
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voltages  of  hydrophones  at  locations  l,  2,  and  3.  respectfully.  The 
hydrophones  are  in  direct  communication  with  the  T8L  and  the  signals  are 
represented  as  a  linear  combination  of  various  source  terms  (Equations  2.1, 
2.2.  and  2.3 ). 

a(r,  .e.t)  *  a^Cr,  .9.0  ♦  aU'(r,  .9.0  ♦  atb)(rj  ,9,0  ♦  ae(r,  .9.0  (2.1) 

b(r,  ,9,t)  =  bAc(r,  .9.0  ♦  bu  (r,  .e.t)  ♦  btb,(r,  .e.t)  ♦  be(r,  .e.t)  (2.2) 

c(r,  ,9,t)  =  cAc(r,  ,9,t)  ♦  cu  (r,  .e.t)  ♦  ctbt(r,  ,9,t)  ♦  ce(r,  .e.t)  (2.3) 


Figure  2.1:  Array  and  Sensor  Locations. 


In  the  above  equations  subscripts  Ac,  u\  tbl,  and  e,  refer  to  individual 
sources  that  describe  the  total  signal  measured  at  I.  2.  and  3  due  to  the 
facility  generated,  vibrationally  induced,  TBL  pressure  fluctuations  and 
electronic  noise,  respectively.  The  outputs  u'a(t).  u'b(t),  and  u'c(t)  describe 
the  time  varying  vibrations  measured  by  accelerometers  located  on  the 
tunnel  wall  near  each  respective  hydrophone  location.  Finally,  all  outputs 
are  considered  to  be  stationary,  random  processes  of  zero  mean. 
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2.3  node!  AssumD.tions_an<l  Considerations 

The  foremost,  of  the  model  assumptions  is  that  all  facility  generated 
noise  propagates  as  plane  waves,  constant  in  amplitude  and  phase  in  any 
plane  perpendicular  to  the  flow  in  the  test  section.  If  all  major  sources  of 
the  acoustic  or,  propagating  pressure  field  are  located  several  diameters  up 
or  downstream  from  the  test  section  then— for  frequencies  below  the 
cutoff,  jco  (Appendix  2)— this  foremost  assumption  should  be  valid.  This 
premise  implies  that  by  simply  subtracting  any  pair  of  hydrophone  outputs 
(e.g.  aCrj.e.t)  -  b(rj,e,t)  ),  virtually  all  acoustic  or.  propagating  pressure 
contributions  can  be  eliminated.  Also,  if  the  primary  objective  is  to 
determine  the  spectrum  of  the  TBL  wall  pressure  fluctuations  then  the 
following  restrictions  on  the  measurement  system  apply: 

1)  All  hydrophones  must  have  the  same  receiving 
sensitivities 

2)  All  hydrophones  must  be  electronically  matched 
in  phase 

3)  Hydrophone  separation  distances  must  be  greater 
than  the  maximum  correlation  length  of  the  TBL 
pressure  fluctuations  (generally  accepted  to  be  on 
the  order  of  the  boundary- 1 3yer  thickness). 

If  either  of  the  first  two  criteria  are  not  true,  then  the  acoustic  or, 
propagating  component  of  the  pressure  field  will  not  be  effectively 
cancelled.  In  addition,  if  item  three  is  not  met  then  an  attenuation  of  the 
local  TBL-induced  pressure  signal  will  result.  Therefore,  all  hydrophones 
are  assumed  matched  in  phase,  equal  in  receiving  sensitivity  and  their 
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separation  distances  are  maximized  to  allow  negligible  attenuation  in  the 
local  TBL  pressure  fluctuation  levels  ( see  Chapter  4  for  further  discussion 
on  correlation  length  scales).  Summarizing  the  above  concepts  suggests 
that  the  following  source  related  terms  are  either  correlated  or, 
uncorrelated: 

1)  Acoustic  source  terms  are  uncorrelated  with  TBL  source 
terms 

2)  TBL  pressure  fluctuations  at  a  given  location  are 
uncorrelated  with  those  at  other  locations 

3)  Acoustic  terms  are  fully  correlated  between  hydrophones 

4)  Electronic  source  terms  are  uncorrelated  between  hydrophones 
and  with  ail  other  non-electronic  terms. 

2.4  AQalusis :  .the  finite  Foucier  Transform 

The  Fourier  transform  of  a  nonperiodic  function  x(t)  can  be  represented 
as  (note  that  the  upper  case  lettering  refers  to  the  transformed  quantity). 


X(f)  *  F{  x(t) 


oo 

) 3 J  x(t)  e"2^^  dt. 
-00 


(2.4) 


It  is  not  possible  to  obtain  data  over  an  infinite  record  length,  therefore 
the  Fourier  transform  of  x(t)  must  be  approximated  by  integrating  over  a 
finite  record  of  length  T.  For  example,  x(t)  is  approximated  as, 

x(t)  =  (  xk(t)  [  H(t)  -  H(t-T)  ] }  0<  k  <  n  (2.5) 


A'Z 

.  * .  -  w >  i 

w.-.y.* 

V-V'C-' 
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— n/vt  jv t*vt  ■»“ 
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where,  nT  is  the  total  record  length  of  x(t),  k  is  an  individual  record  of 
length  T,  and  H(t)  is  the  Heaviside  Function  defined  as, 


H(t)  =  0.  for  t<0 

H(t)=  1.  for  t  >  0. 

Substituting  Equation  2.5  into  Equation  2.4  yields  the  finite  Fourier 
transform  of  the  klh  record  of  xj<(t)  ( -  x(t,k) ). 

00 

X(f,  T,  k)  »  F(  x(t,  k) )  -  j  x(t,  k)  (  H(t)  -  H(t  -  T)  ]  dt 

-00 

00  00 

•  |  xU.Ue-^dt-  {  xa.kJe-^JKdt 

0  T 

T 

■  J  x(t,  k)  e'217^  dt.  (2.6) 

0 

2.5  IBL  Spedimfrem  Cross-Difference  Signals 
The  finite  Fourier  transform  of  the  difference  signals  between  two 
pairs  of  hydrophones  is  the  autospectral  density  of  the  TBL  pressure 
fluctuations.  This  analysis  is  shown  as  follows. 

From  Equation  2.6,  the  finite  Fourier  transform  of  the  difference  signals 
between  the  hydrophone  outputs  at  locations  1  and  2.  and,  I  and  3  can  be 
constructed  and,  is  represented  by  Equation  2.7.  However,  in  Equation  2.7 
the  dependency  on  the  spatial  coordinates  r,  and  9  have  been  suppressed  for 
convenience  (i.e.  the  Fourier  transform  is  evaluated  where  time  is  the 
variable). 
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[  A(f,  T,  k)  -  B(f .  T,  k)  ] 


T 

{  (taAc<t)  ♦  aU'(t)  ♦  atb|(t)  ♦  ae(t)l 

-lbAc(t)  ♦  V(t>  btb,CO  ♦  be(t)l)e'27TJ^  dt. 


and, 

l  A(j,  T,  k)  -  C(j,  T,  k) 


(2.7) 

T 

\  ♦  a^(t)  ♦  att),  (t)  ♦  ae(t)l 

“(cA<Xt)  +Cy*(t)  ♦  Cg(t))}e"27rilt  dt. 


(2.8) 

In  terms  of  their  individual  transformed  quantities  Equations  2.7  and  2.8 
become,  where  the  arguments  have  been  suppressed  for  simplification. 


DA-B(f-  T,  k)  S  l  A(f,  T,  k)  -  B(f,  T,  k)  1  •  Atb1  -  Btbl  -  AAc  -  BAc 


and, 


*  Ay'  -  BU' ♦  Ae  -  Be  (2.7a) 


DA-CCJ,  t*  k)  *  l  A^»  T' J  ■  Atbl  "  ctbl *  aAc  *  cAc 

♦  Ay'  ”  Cy'  *  Ag  Ce.  (2.8a) 


Appendix  1  shows  that  the  autospectral  density  function.  Gxx(j)  and 
cross-spectral  density  function  GXy(f) can  b 2  represented  by  the  respective 


products  of  their  Fourier  transformed  time  functions  x(t)  and  y(t). 
Therefore,  the  cross-spectral  density  function  between  the  Fourier 
transformed  difference  signals  (Equations  2.7a  and  2.8b)  is  given  by, 

GA-B,A-dl>  “  2  llm  £((A  -  B)*  (A  -  C)]/  T.  (2.9) 

T 

Expanding  the  terms  inside  the  expectation  operation  of  Equation  2.9  gives. 

Ga-B,A-c(I>  5  lim  I  E  < lA-,0,  -  B-foi  *  A*u-  -  B»u-  •  A*e  -  B*cl 

T~»00  _ 

T  (A^|  -  Ctb,  >  AU'  -  CU'  +  Ae  -  Ce]  ). 

(2.10) 

In  F  n  2.10  the  expected  value  of  all  uncorrelated  terms  will  tend  to 
zero  so  that, 

GA-B.A-C<f)  s  2  lim  I  {E{A*tb,  Atbl]  -  E[AV  AU'l  -  E(A*V  Cy'l  ♦  E(BV  Cy-1 

T-»00  - 

T  -ElB V  AU1  -  ElB V  V  1  *  EfAwu'  Atb|  J  -  E[A  V  Ctb)  I 

♦E(BVctb|]  “  E(B*V  Atb|]  ♦  E[A**(b|  A^]  -  E[A**tb|  Cu] 
♦ElB-tb|  Cy'l  -  E[B**|b,  Ay']  ♦  E(A**g  Ag] ).  (2.11) 

Solving  Equation  2.11  for  the  expected  value  of  the  turbulent  boundary 
layer  pressure  fluctuations  at  location  I.  and  by  combining  similar 
quantities  it  easily  shown  that, 
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Ga^iO)  =  6A-8,A-<$)  -  2  lim  E((  AU'  -  By' )"(  AU'  -  CyO)  /  T 

T-»oo 


-  2  llm  E[(AU-  -  Bu-)*(Albl  -  Ctb))l  /  T 

J  -*00 


-  2  lim  Et(  Atb,  -Btt)|)*(Ay'  -  Cy')]  /  T  -  Ga  (J). 

T-*eo 

(2.12) 


2.6  Analysis  of  Vihrationatly  Induced  Terms 

Suppose  that  a  hydrophone's  output.  y(t)  is  due  only  to  its  local 
acceleration.  u'(t).  Then  with  all  extraneous  noise  absent,  the  hydrophone's 
output  can  be  modeled  as  an  ‘‘ideal.''  single  input/output  system,  as 
depicted  in  Figure  2.2. 


u'(t),U'(f) 


h(t),  H(|) 


»  y(t),  Y(f) 


Figure  2.2:  Ideal  Single  Input/Single  Output  System. 

For  this  ideal  single  input/output  system  y(t)  is  given  by  the 
convolution  of  the  impulse  response  function.  h(r)  with  the  system  input 
u'(t)  or,  using  the  Fourier  transform, 


Y(j)  =  H(j)U'(j). 


(2.13) 
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where  the  upper  case  lettering  refers  to  the  usual  Fourier  transformed 
quantities.  From  Equation  2.13  it  follows  that. 


Y"(j)  =  H*(j)  U'«(j) 

Y“(  f)  Y(f)  =  H*(f)  H({)  U'"(f )  U'(f)  (2.14) 

and, 

U'*(f)  Y(j)  =  H({)  U'*({)  U'(f).  (2.15) 


Multiplying  Equations  2.14  and  2.15  by  2/T,  letting  T  tend  to  infinity,  and 
taking  the  expected  values  give  the  usual  spectral  relationships. 


Gyu(n  =  1  H(f)l2euu.(f) 

(2.16) 

Gu'y(f)  s  H(f)  GU'U'(f). 

(2.17) 

Clearly,  H(j)  is  a  complex  function  and  for  convenience  will  be 
represented  as, 


H(f)  =  |  H(j)  |2e-j©CP 

where, 

e(|)  =  tan"1  [  lm(  Gxy(()  /  Ret  Gxy(f) )  ]. 


Using  these  and  previous  results,  the  transfer  functions  between  the  locally 
measured  accelerations  (u'a(t),  u'j>(t),  u'^t))  and  the  respective 
hydrophone's  vibrational  outputs  can  be  constructed  and  are  given  by  the 
following  spectral  relationships, 


HA(j)  =  Ay'(J)/  U'A(f) 

(2.18a) 

HB(f)  =  Bu'(f)/  U'B(f) 

(2.18b) 

HC(|)  =  Cu(l)/  U  cd). 

(2.18c) 

In  practice,  these  transfer  functions  are 

measured  under  a  no-flow 

condition  (all  non-vibrational  source  terms  are  forced  to  zero)  using  some 
form  of  mechanical  excitatioa  Solving  Equations  2.18  (where  frequency 
variables  have  been  suppressed)  for  their  vibrational  outputs  result  in, 

V  (D=  u-  ha 

(2.19a) 

Bu-<|)  =  U-  HB 

(2.19b) 

Cu-(f)  =  IT  HC 

(2.19c) 

Substituting  Equations  2.19  into  Equation  2.12  yields, 

6Atbl^  s6a-b,a-cW 

-  2  lim  E((  U'A  Ha  -  U'B  Hb)"  (  U'A  HA  -  U'c  Hc)  ]  /  T 

T  -*oo 

-2  lim  £((  U' A  Ha  -  U'b  HB)»  (Att>|  -  Ctt), )  1  /T 
T  ”»PO 

-  2  Mm  El  (Atb|  -B,„| )«  ( ll'A  Ha  -  U'c  HC)  ]  /  T  -  G.  (J). 

T-.00 


In  tne  above  equation,  terms  that  involve  the  turbulent  boundary  layer 
pressure  fluctuations  are  not  directly  measurable;  however,  there  exists 
quantities  that  involve  these  source  terms  which  can  be  measured.  By 
substituting, 


A<b|  (f.T.k)  ”  Cjjjj  (f.T.k)  =  A  “  C  -  (Ay*  ~  CyO  _  (Ag~  Cg) 
and, 

Atbl(f,T.k)  -  Btb|  (f.T.k)  =  A  -  B  -  (Ay'  -  BU')  -  (Ae  -  Bg). 


into  Equation  2.21,  we  find, 


GAtb|0) s  6A-b,A-C<1) 

♦  2  lim  1  El(U'AHA-U'BHB)-(U'AHA-U'cHc)l 
T  -*oo 

T 

-  2  lim  1  E[  (  U'A  Ha  -  U'B  HB)*  (A  -  0  ] 

T  -*oo 

T 

-  2  lim  1  El  (A  -  B)"  (  U'A  HA  -  U'c  Hc)  ]  "  ^P. 

J  -*oo 

T  (2.22) 

If  the  measured  complex  transfer  functions  in  the  above  equation  are 
not  equal  in  magnitude,  and  phase  over  the  desired  frequency  range  (e.g. 
Ha*Hb*Hc)  then  al!  expected  value  terms  in  Equation  2.22  must  be  expanded 

to  yield  the  desired  measurable,  spectral  quantities.  Expanding  Equation 
2.22,  and  taking  the  limit  of  individual  terms  as  T  tends  to  infinity  will 
then  yield  the  following  expression  for  the  TBL  pressure  fluctuations  in 
terms  of  other  measurable  spectral  quantities, 


6Atb,(|)  =  GA-B,A-C  "  H"AGU'AA  4  H"agU'aC 

+h*b  gitba  •  h*b  gubc  "  ha  gau' a 

4HC  gAJ'C  4  hA  gBU'a  “  HCgBU'c 

4Iha|2(5U  aU'a  '  H*AHCGU'AU'c 

■  h"b  hA  gU'BU' a  4  H*B  hC  gU'BU'c  *  6v  <2-23> 

However,  if  the  complex  transfer  functions  are  equal  in  magnitude  and 
phase  (Ha=Hb=HcsH)  for  all  frequencies  of  interest  then.  Equation  2.22 
reduces  to. 

GAtb|W  =  ga-B.a-C  * °a2  ' H"  %A-UB.A-C 

"HGA-B.U'a-U'c  *  lHf  %a“U'b- UA"U  C'  (2'24) 

2.7  Estimation  of  viDratlonal  Spectral  Quantities 

It  is  possible  to  gain  insight  into  the  relative  magnitude  of  the 
vibrationaily  related  terms  in  Equation  2.24  without  actually  measuring 
them.  Consider  the  system  model  shown  in  Figure  2.3. 

In  this  model  y(t)  represents  the  output  due  to  the  difference  between 
two  hydrophone  signals,  and  n(t)  is  that  part  of  y(t)  due  to  the  difference 
between  TBl  source  terms  only.  Finally.  x(t)  represents  the  signal  due  the 
difference  between  the  respective  hydrophone's  accelerometers  outputs. 

The  true,  uncontaminated  vibrational  difference  signal  output.  v(t)  of 
the  two  hydrophones  is  therefore  obtained  by  the  convolution  of  h(r)  with 
x(t).  All  electronic  source  terms  are  assumed  negligible  in  this  model. 
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x(t)  =  u'a- 


y(t)  =  a  -  b 


Figure  2.3 :  Vibrational  Noise  Model 
From  Figure  2.3  It  Is  clear  that, 


Gyy  (f)  =  GnnW  *  Gy\/f) 

Gy\<f)  s  |  H(f )  | 2  GxxCP 

Gxy(f)  =  H(f)  Gxxd>- 

The  ordinary  coherence  function  between  x(t)  and  y(t)  Is  defined  as, 


2  «GXY(^I 

4(f)  =  -  <  1. 

Gxxd) 


From  the  above  relationships  It  follows  directly  that, 


Gyv(f)  = 


GXY(^I2  .  . 2 


G2xxd) 


GxxW  =  ^YY^)- 


The  above  equation  Is  termed  the  coherent  output  power  (COP)  spectrum 
between  y(t)  and  x(t)  and,  represents  the  power  spectrum  of  v(t)  with  the 
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linear  effects  of  n(t)  removed.  Therefore,  if  the  COP  spectrum  of  y(t)  is 
much  less  than  the  autospectrum  of  y(t)  (6  yy  (f)).  the  dominant  source  term 
is  n(t)  or.  that  part  of  y(t)  due  to  the  TBL  pressure  fluctuations.  In  other 
words,  if  the  ordinary  coherence  function  between  x(t)  and  y(t)  is  much  less 
than  unity,  the  output  y(t)  is  dominated  by  the  local  TBL  pressure 
fluctuations  and  all  vibrational  terms  in  Equations  2.24  and  2.25  may  be 
neglected.  If  this  is  true,  then  Equation  2.24  reduces  to. 

%|(1)  =  0A-B.A-C(|)-S(f)-  (Z25) 

This  equation  for  Gtbid)  represents  the  true  turbulent  boundary  layer  wall 
pressure  spectrum  with  all  extraneous  noise  removed.  The  electronic  noise 
spectrum.  is  easily  determined  by  a  direct  measurement  of  the  output 
at  sensor  a.  with  the  facility  off  and  is  usually  negligible. 

2.8  Cross-Spectral  Measurements:  Acoustic  and  Vibratory  Noise 

Measuring  the  cross-spectral  density  function  between  any  two 
hydrophones  will  yield  only  those  noise  sources  that  are  correlated 
between  the  two  points  of  measurement.  In  other  words,  the  local  TBL 
pressure  fluctuations  will  be  conditioned  out  of  the  cross-spectrum 
assuming  the  measurement  points  are  not  In  close  proximity.  For  example, 
consider  the  cross-spectrum  between  hydrophone  signals  a  and  b, 


6abW*2,,it1  1  El (Atb!  *  aAc *  AU' 4  Ae>* (®tbl  *  bAc 4  BU'  * 
T  “♦00  — 

T 


(2.27) 


Expanding  the  expectation  operation  and  keeping  only  non-zero  terms  gives, 


Gab<J)-  2  iim  1  3(j']  ♦  ElA^'B^jl  ♦  E(A*y'  Bu'l 

J  -»oo - 

T  4  ElA^AcBAc]  ♦  E(A«Ac  B(jO  ♦  E[A*U'BAC] }. 

(2.28) 

Recalling  that. 

Att)|(f)  °  A  -  Ay'  -  Ae  -  AAc 
Btb|(|)  -  B  -  By-  -  Be  -  BA. 

Substituting  these  expressions  Into  Equation  2.28  and  collecting  like  terms 
yields  the  following  expression, 

G^p^llm  I  {  E(A*B^|  ♦  ElA^B]  -  ElA*^'1  *  E<A*AcBAc  >* 

T-*oo  — 

T  (2.28) 
From  Equation  2.19  the  above  equation  reduces  to. 

=  hB  gAU'q  *  A  eLT AB  ’  H“ A  hB  GU  AU'B  *  °*A&AC  (2‘30) 

The  fourth  term  on  the  right  hand  side  (R.H.S.)  of  Equation  2.30 
represents  the  true  propagating  acoustic  pressure  field  and  propagates  as 
plane  waves  throughout  the  test  section  (i.e.  constant  in  magnitude  and 
phase).  All  other  terms  on  the  R.H.S.  of  Equation  2.30  represent  a 
pseudo-acoustic  field  due  to  the  local  wall  and/or  transducer  vibrations. 


?g  Fstimation  of  Vibrational  Terms 

Again,  using  an  analysis  that  is  similar  to  the  one  used  previously,  it  is 
possible  to  determine  the  relative  magnitude  of  the  vibrational 
contamination  present  in  the  cross-spectrum  of  sensors,  a  and  b  (without 
measuring  the  transfer  functions  HA(j)  and  Hg(j)).  Consider  the  model 
shown  in  the  Figure  2.4  . 


n(t)  =  atb1(t) 


y(t)  =  a(t) 


Figure  2.4:  vibrational  Noise  Model. 


In  Figure  2.4,  y(t)  represents  the  output  from  a  single,  hydrophone.  n(t) 
is  part  of  y't)  due  to  the  local  TBt  pressure  fluctuations,  v(t)  is  the 
hydrophone's  vibrationally  induced  output  ano  u'a(t)  is  the  corresponding 
accelerometer’s  vibrational  output.  The  electronic  noise  source  term  is 
neglected.  Measuring  the  coherent  output  power  (COP)  spectrum  between 
y(t)  and  u'(t)  will  then  determine  the  amount  of  vibrational  contamination 
present  in  y(t).  If  the  COP  spectrum  is  much  less  (<lOdB)  than  the 
autospectrum  of  y(t),  all  vibrational  terms  in  Equation  2.30  can  be 
neglected.  If  this  is  true  then  Equation  2.31  will  reduce  to. 


=  gAAcBAc  =  G  Acoustic 


The  above  equation  represents  all  facility  generated  noise  that  propagates 
as  plane  acoustic  waves  in  the  test  section. 
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Chapter  3 

PRESSURE  SENSORS  AND  CALIBRATION  TECHNIQUES 
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3.1  introduction 

In  Chapter  2  a  measurement  technique  was  developed  that  describes  how 
an  axisymmetric,  pressure  sensor  array  can  discriminate  between  the 
propagating  (acoustic)  and  nonpropagating  (turbulent  boundary  layer.  TBL) 
pressure  fluctuations  in  a  fully  developed,  turbulent  pipe  flow.  However,  a 
prerequisite  and  critical  aspect  of  this  measurement  technique  is  that  1) 
the  phase  difference  between  pressure  sensors,  and  2)  the  receiving 
sensitivities  of  each  sensor  match.  Moreover,  the  facility  in  which  the  TBL 
pressure  measurements  are  to  be  made  uses  glycerine  as  the  working  fluid. 
Therefore,  possible  differences  between  water  and  glycerine  hydrophone 
sensitivities  may  exist  and  is  investigated. 

This  chapter  describes  the  type  of  pressure  sensors  used  and  also,  the 
calibration  schemes  employed  in  clarifying  the  above  sensor  requirements. 

3.2  Pressure  Sensor  Description 

In  order  to  minimize  possible  phase  and  sensitivity  differences  between 
pressure  sensors,  high  quality,  resistive-type  sensors  from  the  ENDEVCO 
Corporation  were  chosen.  Resistive-type  sensors  have  a  purely  real  input 
impedance  and  should  theoretically  introduce  no  phase  delays  in  the  output 
signal.  Also,  high  quality  pressure  sensors  should  have  a  flat  frequency 
response  (to  first  resonance)  and  more  importantly,  the  sensitivities 
should  vary  little  from  sensor  to  sensor.  Another  aspect  of  the  pressure 
sensors  that  was  mentioned  in  Chapter  1  was  that  the  sensor's  facial 
diameter  has  to  be  small,  in  relation  to  the  smallest  size  of  the  TBL 


.*>>1  "  v.'<i 


pressure  producing  eddies.  If  the  sensor's  diameter  is  large,  then  local, 
spatial  averaging  of  the  pressure  producing  eddies  will  result  producing  a 
mean  output  voltage  that  approaches  zero  (see  Figure  3.1). 


j  Figure  3.1:  TBL  Pressure  Fluctuations  and  Sensor  Size. 
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)  Therefore,  suDminiature  pressure  sensors  with  a  casing  diameter  of  1.57 

millimeters  (mm)  were  selected  (ENDEVCO,  Model  Number  8514-10).  It 
should  be  noted  however,  that  the  active  sensing  elements  of  each 
transducer  are  much  smaller.  Galib  and  Zandina  (1984)  have  used  these 
types  of  transducers  for  TBL  wall  pressure  measurements  and  have 
concluded  (through  testing)  that  the  active  area  is  probably  as  small  as  0.5 
mm.  All  sensors  are  rated  to  6.8944*104  N/m^with  a  corresponding  full 
scale  output  of  300  mV.  Also,  each  sensor  can  withstand  transient  inputs 
of  6.8944*10^  N/m2  without  damage  and  outputs  remain  linear  for 
pressures  up  to  three  (3)  times  the  full  scale  (F.S.)  output.  Due  to  the  harsh 
thermal  environment  that  exists  in  the  Glycerine  Tunnel,  each  transducer 


was  thermally  compensated  for  temperatures  between  -17.7®  and  93.3°C. 
Each  sensor  was  specially  ordered  with  a  double  coating  of  “Parylene  C.‘  (a 
factory  sealant)  that  provided  more  than  adequate  sealing.  The  transducers 
consist  of  a  zero  balanced,  four -active-arm,  strain  gage  bridge,  diffused 
into  a  molded  silicon  diaphragm  (see  Figure  3.2).  Thermal  compensation 
resistors  (Rj,  R2)  and  the  complete  bridge  circuit  are  all  housed  in  the 
same  nickel-iron  alloy  case.  Powering  the  transducers  with  a  10  volt  DC 
source  is  recommended;  however,  appreciable  shift  in  output  sensitivity  is 
not  observed  at  even  half  the  recommended  supply  voltage  (  see  Figure  3.3). 


Figure  3.2:  Schematic  Diagram  of  the  4-Arm  Bridge  of  the  Sensors. 


The  power  supply  requirement  was  fulfilled  with  eight  1. 25-volt, 
dry-cell,  rechargeable  nickel-cadmium  batteries  (10-volt  total  supply).  A 
major  advantage  of  batteries  over  commercial  DC  power  supplies  is  the 


Figure  3.3:  Typical  Transducer  Sensitivity  Shift  (mV/Pa/Input  Volt) 

Due  to  Varying  Supply  Voltage. 

elimination  of  all  60  cycle  related  noise.  The  current  needed  to  power  each 
sensor  was  small  ( 6mA/sensor),  enabling  the  powering  of  all  three  sensors 
simultaneously.  The  power  supply  unit  W3S  tested  and  found  capable  of 
supplying  the  required  input  voltage  for  up  to  seven  (7)  hours  of  continuous 
operation  without  recharging  (Figure  3.4).  Factory  supplied  free-field 
pressure  sensitivities  for  each  sensor  and  other  typical  performance 
specifications  are  given  below  in  Tables  3.1  and  3.2  (note  that  all  factory 
specifications  are  with  a  10-volt  supply).  Shown  in  Figure  3.5  are  the 
physical  dimensions  of  a  typical,  8514-10  unit.  A  removable  screen  that 
consists  of  eleven  extremely  small  slots  protects  the  silicon  diagram  from 
undesired  impacts  and  is  not  shown.  However,  the  transducer  properties 
and  performance  are  uneffected  by  the  screen’s  presence. 
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Table  3.2:  Endevco  Typical  Specifications.  Series  8514  Transducer 
with  a  10-Volt  Power  Supply. 


Performance _ 

Range 

Resonance  Frequency 
linearity  @  FSO 
Linearity  @  3  x  FSO 
Hysteresis 

Zero  Shift  @  3  x  overrange 
Zero  Shift  w/  temperature 
ref.  23.9  *  C 

Sensitivity  Shift  w/  temp. 

ref.  23.9 #  C 
Warm  up  Time 
Acceleration  Sensitivity: 
lateral 
Longitudinal 
Burst  Pressure: 

Burst  Differential 


85M--1Q  .  .  . 

0  -  6.8944  *  104  N/m2 
140,000  Hertz 

♦  0.5  %  of  FSO 

♦  1.0  %  of  FSO 
0.1  %  of  FSO 
0.1  % 

t  3  %  of  FSO  @  max.  rated 
temp. 

♦  4  %  FSO  9  max.  rated  temp. 

15  seconds 

0.4137  N/m2  /g 
10342  N/m2  /g 

♦  6.8944  x  IQS  N/m2 


Figure  3.5:  Physical  Dimensions  of  Subminiature  Transducers. 


30 


3.3  Calibration  .and  Phase  Measurement  Techniques 
Considerable  attention  was  given  to  the  various  methods  of  transducer 
calibration.  However,  for  reasons  of  simplicity,  reliability  and  relative 
accuracy,  the  'comparison  calibration,'  method  that  utilized  a  known 
standard  transducer  was  employed.  The  sensitivity  of  an  unknown 
transducer  is  found  by  subjecting  both  the  unknown  and  the  standard  to  the 
same  incident  pressure  field  and  recording  the  output  voltages  "of  each  unit 
at  various  frequencies  (j).  Then,  the  sensitivity  of  the  unknown 
transducer  is  computed  from, 


vunkr>ownW 

^unknownW  s  MstandardW  '  ^  ^  ^ 

standard*^ 

Phase  calibration  imposed  an  additional  restriction  on  the  calibration 
device.  In  order  to  keep  the  phase  of  the  incident  pressure  wave  at  the 
face  of  each  sensor  the  same,  a  plane  wave  sound  field  is  needed.  This  was 
accomplished  by  utilizing  a  thick-walled,  glass  tube  as  the  calibrator.  For 
a  cylindrical,  rigid-walled  tube,  the  cutoff  frequency  for  plane,  axially 
symmetric  wave  propagation  can  be  obtained  and  is  shown  in  Appendix  2 
to  be, 


jco- 0.293  C/R.  (3.2) 

In  the  above  equation,  c  represents  the  sound  speed  of  the  calibration 
medium  and  R,  the  radius  of  the  tube.  In  the  case  of  the  5.08  centimeter 
(I.D.)  tube  used  for  these  calibrations  a  corresponding  cutoff  frequency 


was  computed.  If  water  filled,  the  cutoff  frequency  of  the  tube  is  17,000 
Hertz  and.  22.800  Hertz  if  glycerine  filled.  These  frequencies  represent  an 
upper  limit  to  which  calibrations  and  moreover,  phase  data  are  accurate. 

3.3.1  Calibration  Tube  Assembly 

The  calibration  tube  and  assemblies  (Figure  3.6)  were  vertically 
positioned  and  a  piezo-ceramic  source  (projector)  placed  and  sealed  into 
the  lower  end  of  the  tube.  A  large,  cylindrical  shaped  transducer  was 
chosen  as  the  source.  Its  facial  diameter  was  4.76  centimeters  (cm)  and 
had  a  length  of  3.81  cm.  Although  the  transmitting  response  of  the  source 
transducer  was  not  critical  because  a  comparison  calibration  method  was 
used,  the  received  signal  from  the  source  was  always  at  least  10  dB  above 
the  background  noise  level  ( i.e.  noise  present  at  receiver  with  source  off). 
The  opposite  end  of  the  tube  was  customized  to  accommodate  either  all 
three  small  area  transducers  or,  the  standard  LC-10  hydrophone. 

The  small  area  transducers  were  all  housed  in  a  single  assembly  unit  and 
were  circumferentially  placed  120-degrees  apart.  This  assembly  was 
placed  in  the  far  end  of  the  calibration  tube  and  rested  on  a  vibration 
isolating  o-ring.  The  LC-10  hydrophone  had  a  similar  assembly  and  could 
also  be  positioned  in  the  far  end  of  the  tube.  Both  assemblies,  when 
installed  were  located  10  feet  from  the  source. 

3.3.2  Water  Versus  Glycerine  Calibration  Mediums 

Comparison  calibration  techniques  reference  all  measurements  to  the 
original  calibration  medium  of  the  standard.  Therefore,  because  the 
standard  (LC-10)  used  in  these  calibrations  was  (factory)  calibrated  in 
water,  the  calibration  results  of  8514-10  transducers  would— regardless 
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of  calibration  medium— be  referenced  to  water.  Since  the  characteristic 
acoustic  impedances  between  water  and  glycerine  (pcwatcr/  pCgiycerine  5 
0.592)  differ,  possible  discrepancies  in  calibrations  involving  these  two 
mediums  were  considered. 

Air  and  water  methods  of  transducer  calibration  will  yield  the  same 
results  at  audio  frequencies  if: 

1)  the  transducer's  radiation  impedance  can  be  neglected  and, 

2)  the  transducer  is  small  enough  compared  to  the  acoustic 
wavelength  of  the  medium  so  that  diffraction  effects 
can  be  neglected. 

The  above  criteria  suggest  that  if  ka«l,  where  k  is  the  acoustic 
wavenumber  and  a,  the  transducer  radius  then,  the  transducer's  sensitivity 
should  be  Independent  of  the  calibration  medium. 

To  test  this,  the  output  voltage  of  a  8514-10  transducer  was  recorded 
(with  a  white  noise  source  Input)  when  the  calibration  tube  when  filled 
with  both,  water  and  then,  glycerine  mediums.  Results  shown  In  the  above 
Figure  3.7,  suggests  no  difference  in  calibration  will  occur  over  the 
specified  frequency  range  if  water,  rather  than  glycerine  were  used  as  the 
calibration  medium.  Therefore,  because  of  its  relative  handling  ease  water 
was  chosen  as  the  calibration  fluid. 
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Figure  3.7:  water  and  Glycerine  Transducer  Output  Voltage  versus 
Frequency. 

3.4  Data.AcqulaltlQn.afxl  Reduction 

Calibration  and  phase  data  were  needed  over  a  broad  range  of 
frequencies  (0  to  10  kHz).  This  was  accomplished  by  using  broad-band 
noise  as  Input  to  the  source.  In  this  way,  a  quick  method  of  calibration 
over  the  desired  frequency  range  was  possible  using  FFT  processing. 

The  output  of  a  General  Radio  random  noise  generator  was  amplified  and 
used  to  drive  the  source  transducer.  This  input  was  monitored  on  a  dual 
channel  Spectral  Dynamics  360  spectral  analyzer  (5D360)  and  was  flat  over 
the  frequency  range  of  interest  (Figure  3.8).  Outputs  from  the  hydrophones 
(small  area  and  LC-10)  were  individually,  bandpassed  filtered  and  amplified 
via  a  Brookdeal.  differential  filter/amplifier.  Half  power  points  of  the 
highpass  and  lowpass  filters  were  set  at  1  and  10,000  Hertz,  respectfully. 
These  outputs,  were  then  processed  on  the  SD360  spectral  analyzer.  All 
analyzer  data  was  averaged  512  times  over  a  frequency  range  of  0  to 
10,000  Hertz. 
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Figure  3.8:  Source  Input  Voltage  Spectrum. 


This  frequency  range  of  calibration  set  the  corresponding  bandwidth  of 
the  Hanning,  analysis  window  to  20  Hertz.  After  processing,  all  data  were 
sent  to  the  ARl/VAX  782  computer  for  reduction.  Figure  3.9  depicts  the 
data  acquisition  system  employed  throughout  the  calibrations. 

Sensitivities  of  the  small  area  transducers  were  obtained  by  recording 
their  output  voltages.  Then,  for  the  same  source  input,  the  small  area 
transducers  (8514-10)  were  replaced  with  the  standard  LC-10  and  its 
output  voltage  was  recorded.  These  data  were  then  processed  on  the  VAX 
computer  according  to  equation  3.1,  yielding  the  desired  sensitivities. 

The  phase  relationships  between  the  small  area  transducers  were 
calculated  directly  by  the  SD360.  That  is,  the  phase  differences  were 
obtained  by  measuring  the  cross-spectral  density  function  (G^g(j))  between 
two  pressure  sensors.  The  phase  difference  was  tr)en  calculated 

from, 
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Figure  3.9:  Data  Acquisition  System  and  Tube  Calibrator. 

Phase  differences  were  then  plotted  on  the  5D360's  analog  plotter  for  all 
three  cross-spectral  combinations  of  the  small  area  transducers. 

3.5  figSilUS. 

The  measured  sensitivities  or  the  small  area  transducers  are  shown  as  a 
function  of  frequency  in  Figures  3.10,  3.11  and  3.12.  Due  to  the  relatively 
flat  frequency  response  in  the  measured  sensitivities  of  all  three  pressure 
sensors,  averaged  values  are  given.  These  averaged  values  of  the  small  area 
transducers  are  all  within  0.5  dB  of  the  factory  supplied  calibrations. 
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Shown  in  Figures  3.13,  3.14  and  3.15  are  the  measured  phase  differences 
between  the  small  area  transducers.  All  plots,  for  the  frequency  range 
specified  show  a  negligible  (phase)  shift  in  the  received  signal  between 
paired  sensor  outputs. 

3.6  CQQfli&LQQS 

The  techniques  employed  in  measuring  the  sensitivities  and  phase 
relationships  of  the  small  area  transducers  provided  an  excellent  way  to 
check  the  factory  specifications  and  theoretical  considerations.  The 
measured  sensitivities  between  sensors  varied  only,  0.5  dB  over  the 
frequency  range  of  calibration  and  phase  differences  between  the  sensor 
outputs  were  insignificant.  Overall,  the  sensors  proved  to  be  quality 
devices  and  meet  all  of  the  requirements  for  the  measurement  of  turbulent 
boundary  layer  wall  pressure  fluctuations. 

A  plane  wave  environment  was  achieved  by  utilizing  a  tube  as  a 
calibration  device.  The  calibrator  proved  that  accurate  free-field 
calibrations  and  phase  measurements  are  possible  over  a  wide  frequency 
range  (0  to  10kHz)  without  the  traditional  use  of  anechoic  chambers. 

Also,  because  of  the  tube's  small  volumetric  size,  transducer 
calibration  is  easily  performed  in  any  number  of  test  mediums  (e.g.  water, 
oil,  air,  etc.).  Variations  in  transducer  sensitivity  due  to  different 
calibration  mediums  can  be  investigated,  and  accounted  for  by  simply 
measuring  and  comparing  the  output  voltage  of  the  transducer  when 
immersed  in  both  mediums.  The  sensitivity  of  the  transducer  in  the  new 
calibration  medium  (medium  2)  is  given  simply  by  the  following 
relationship. 


vmedium2  (P 

ttmedium2(P =  nmedium,(P  *  * 

vmediumj  (P 

Although  no  differences  in  sensitivity  between  water  and  glycerine 
calibration  mediums  were  observed  for  the  small  area  transducers  used  in 
this  investigation,  larger  transducers  may  produce  different  results 
because  of  possible  changes  in  impedance  characteristics. 


Figure  3.10:  Measured  Sensitivity  of  Sensor  KW45.  Average  =  -227.6  dt>. 


Chapter  4 

THE  GLYCERINE  TUNNEL  AND  WALL  PRESSURE  MEASUREMENTS 

4.1  introduction 

This  chapter  describes  the  facility  in  which  the  turbulent  boundary  layer 
(TBL)  wall  pressure  fluctuation  measurements  were  made.  These  include 
measurements  of  key  parameters  that  help  characterize  the  TBL  wall 
pressure  spectrum.  This  involved  measuring  the  mean  flow  velocity  (u)  and 
wall  shear  stress  (tw)  as  functions  of  glycerine  temperature.  Also, 
because  empirical  formulae  exist  for  calculating  the  wall  shear  stress, 
theoretical  and  measured  values  of  this  parameter  are  compared. 

This  chapter  also  describes  how  the  TBL  wall  pressure  spectrum  is 
influenced  by  the  Reynolds  Number  of  the  flow.  Reynolds  Number,  is  a 
parameter  that  characterizes  all  types  of  viscous  flows  (internal  and 
external)  and  is  a  ratio  of  a  fluid's  inertial  to  viscous  forces.  Turbulent 
flow  at  the  centerline  of  a  pipe  is  predominantly  influenced  by  the  fluid's 
inertia.  Moving  radially  outward  from  the  centerline,  viscous  forces 
become  more  influential  until,  at  the  wall  the  flow  velocity  becomes  zero. 
The  Glycerine  tunnel  is  a  unique  facility  in  that  the  Reynolds  Number  (Re0  = 
UD/u)  is  varied  by  changing  the  fluid's  viscosity  (o)  rather  than  the  mean 
flow  velocity.  The  viscosity  of  the  glycerine  is  very  temperature  sensitive 
(see  Figure  4.1)  and  therefore,  controlled  by  the  fluid's  temperature. 
Because  of  this,  the  TBL  wall  pressure  spectrum  can  be  measured  over  a 
wide  range  of  Reynolds  numbers  while  maintaining  an  almost  constant  mean 
flow  velocity.  This  provides  a  novel  means  for  investigating  the  scaling  of 
the  TBL  wall  pressure  spectra  in  terms  of  inner  (viscous)  variables.  For 
example,  the  viscosity  (instead  of  the  mean  flow  velocity)  in  this  case 
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represents  the  independent  variable  which  in  turn  affects  rw,  u^,  and  the 
normalized  viscous  time  scale  \ 4  (-  |u/ur). 

Verification  of  one  of  the  model's  foremost  assumptions  namely,  that 
the  TBL  pressure  fluctuations  are  uncorrelated  between  sensors  is  shown. 
This  vastly  improves  the  data  integrity,  especially  at  very  low  frequencies 
(<  10  Hertz),  and  indirectly  proves  that  the  other  TBL  measurement 
assumptions,  described  in  Chapter  2,  are  correct. 

Measurement  of  the  acoustic  background  noise  in  the  facility  was 
investigated  as  part  of  an  assessment  of  the  adequacy  of  the  glycerine 
tunnel  for  future  low- wavenumber,  wall  pressure  spectra  research.  This 
noise  is  shown  to  be  very  intense  for  the  facility  in  its  current 
configuratioa  However,  the  extreme  power  of  the  measurement  technique 
in  its  ability  to  extract  the  desired  TBL  pressure  spectrum  in  the  presence 
of  this  severe,  contaminating  acoustic  noise,  was  demonstrated.  Causes  of 
the  major  noise  sources  is  considered,  along  with  their  possible  solutions. 

Finally,  a  few  suggestions  of  possible  new  areas  of  research  are 
proposed. 

4.2  The  Fxperimental  Facititu 

The  experimental  facility  (Figure  4.2)  is  a  closed  circuit  tunnel  that 
uses  glycerine  as  the  working  fluid.  The  glycerine  employed  in  the  facility 
is  commercially  available  and  in  this  investigation  tested  to  96.5&  purity. 

The  test  section  is  7.6  meters  (24.9  feet)  long  and  has  a  diameter  of 
284.5  millimeters  (11.2  inches).  The  test  section  is  also  honed  to  a 
16-rms  microinch  finish  so  that  all  pertinent  smooth  pipe  formulae  apply. 
Fluid  (glycerine)  is  circulated  by  a  75  kilowatt  constant  speed  (900  rpm), 
centrifugal  pump,  located  in  the  downstream  leg  of  the  test  facility. 
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The  laminar/turbulent  transition  location  of  the  boundary  layer  is  fixed  at 
the  entrance  to  the  test  section  by  a  trip  ring.  Fluid  exiting  the  nozzle 
section  of  the  facility  encounters  an  abrupt  increase  in  cross-sectional 
area  thus,  tripping  the  flow  into  turbulence  at  the  test  section  entrance. 
The  onset  of  turbulence  and  the  location  of  the  transition  zone  were 
investigated  experimentally  by  Bakewell  (1965).  Bakewell's  results 
indicated  that,  for  pipe  diameter  Reynolds  numbers  above  8000,  fully 
developed  turbulent  flow  exists  throughout  the  length  of  the  test  section. 

Reynolds  numbers  are  varied  by  controlling  the  temperature  of  the 
glycerine.  Cooling  is  accomplished  through  a  large,  reverse  flow,  heat 
exchanger  located  in  a  separate  leg,  just  downstream  of  the  test  section 
Heated  fluid  is  discharged  into  the  heat  exchanger  and  after  cooling, 
recirculated  into  the  main  flow  of  the  facility.  Constant  flow 
temperatures  can  then  be  maintained  by  systematically  controlling  the 
cooling  Y£al£L  flow  rate  in  the  heat  exchanger.  Tunnel  temperatures  are 
monitored  by  a  digital  thermometer,  that  is  located  in  the  diffuser  section 
of  the  tunnel. 

The  tunnel  is  covered  with  a  single,  4-inch  wrap  of  fiberglass  insulation 
(tunnel  test  section  has  8  inches).  Bakewell  investigated  the  thermal 
gradient  that  exists  normal  to  the  flow  in  the  test  section.  He  found  that 
the  temperature  gradient  was  minimized  with  the  above  specified 
insulation  thicknesses.  Therefore,  a  uniform  temperature  profile  is 
assumed  throughout  the  tunnel  test  section. 

Because  glycerine  is  an  extremely  viscous  fluid,  small  air  bubbles  tend 
to  become  suspended  in  the  fluid  during  the  tunnel  filling  operation.  To 
alleviate  this  problem,  consecutive  cycling  of  the  tunnel  over  a  period  of 
three  (3)  days  was  required.  The  subsequent  drop  in  the  viscosity  of  the 
fluid  with  heating  allows  all  entrapped  air  bubbles  to  escape  through  the 
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low,  and  high  pressure  observation  domes  of  the  facility.  This  is  then 
verified  by  visually  observing  the  glycerine  (during  tunnel  operation) 
through  the  translucent  domes. 


4.2.1  lyDoeLQgaEiVgiocity 

The  mean  flow  velocity  (U)  was  obtained  as  a  function  of  glycerine 
temperature  by  measuring  the  static  pressure  drop  (PJ-P2)  across  the 
nozzle  section  of  the  tunnel  (Figure  4.3)  for  various  glycerine 
temperatures.  The  Glycerine  tunnel  is  configured  with  small,  static 
pressure  taps  that  begin  at  the  settling  section  and  proceed  in  0.3048 
meter  increments  along  the  entire  length  of  the  test  section.  Two,  small 
lengths  of  capillary  tubing  were  attached  to  two  of  these  taps  (P)  and  P2). 
and  the  respective  high  (Pj)  and  low  pressure  (P2)  sides  of  a  15  PSI.  Bell 
and  Howell  differential  transducer.  The  transducer  was  calibrated  as  a 


Figure  4.3;  Static  Pressure  Measurement  Locations  of  Mean  Flow  Velocity, 


with  Nozzle  Contraction  Ratio  of  16:1. 
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function  of  output  voltage  versus  constant,  input  pressure  differentials  and 
is  shown  below  in  Figure  4.4. 

The  mean  flow  velocity  was  then  determined  by  combining  Bernoulli's 
equation  with  the  conservation  of  mass  law.  This  yielded  the  following 
relationship  between  the  square  of  the  average  velocity  and  the  pressure 
differential  (P1-P2).  In  the  following  equation,  A)  and  A2  are  the 
respective  cross-sectional  areas  at  pressure  locations  P|  and  P2  (see 
previous  Figure  4.3), 


2  2  2  ( P|  -  P2) 

u2  =  U  =  _ t _ 

p  (  I  -  (A2/ A)  )2  1 


(4.1) 


Figure  4.4:  15-PSI,  Differential  Pressure  Transducer  Calibration. 

Due  to  the  large  contraction  ratio  (16.1)  of  the  tunnel  nozzle,  the  square 
of  the  area  ratio  (1:256)  in  the  denominator  of  Equation  4.1  can  be  ignored. 


From  this,  Equation  4.1  reduces  to, 
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The  calculated  mean  flow  velocity  versus  glycerine  temperature  is  plotted 
below  (Figure  4.5)  for  glycerine  temperatures  ranging  from  35°C  to  46.i°C. 


Figure  4.5:  Tunnel  Velocity  Versus  Glycerine  Temperature. 

4.2.2  Wall-Shear  stress  Measurements 
The  wall  shear  stress  was  measured  at  the  same  glycerine  temperature 
(Reynolds  number)  as  that  used  by  Bakewell  and  also,  at  a  number  of  other 
temperatures.  This  was  accomplished  by  measuring  the  static  pressure 
drop  that  occurred  between  two  points  along  the  wall  of  the  test  section. 
The  wall  shear  stress  is  then  obtained  from  the  following  relationship 
(Schlichting  (1979)), 


APD 
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In  Equation  4.3,  D  is  the  pipe  diameter,  and  AP  the  pressure  difference 
measured  between  two  points  separated  by  a  distance,  L  along  the  pipe 
wall.  The  static  pressure  drop  was  obtained  from  two  pressure  taps 
located  5.7,  and  7.3  meters  downstream  from  the  no2Zle  entrance. 
Flexible,  hypodermic  tubing  was  attached  to  the  high  and  low  pressure 
fittings  of  a  5-PS1,  Bell  and  Howell  differential  pressure  transducer  and  to 
the  respective  tunnel,  pressure  taps.  The  pressure  differential  (PpP2) 
was  then  recorded  while  varying  the  tunnel  temperature.  Values  of  the 
measured  wall  shear  stress  were  then  calculated  from  Equation  4.3  and  are 
compared  to  the  theoretical  wall  shear  stress  values  in  Figure  4.6. 
Theoretical  values  were  obtained  from  the  previously  measured  average 
flow  velocities  and  the  following  relationship  for  tW3u  (Schlichting), 

tw  «  0.03955  p  J/4  \)1/4  D  !/4 

Both  theoretical  and  measured  values  are  shown  as  a  function  of  glycerine 
temperature  and  are  in  reasonable  agreement. 


Figure  4.6;  Wall  Shear  Stress  Versus  Glycerine  Temperature. 


Calibration  of  the  Bell  and  Howell  5-PSI  transducer  was  also  done  and  is 
shown  as  a  function  of  millivolt  output,  for  constant  pressure  differential 
input  in  Figure  4.7.  As  with  the  15-PSl  transducer,  an  excellent  linear 
relationship  between  output  voltage  and  pressure  differential  input  was 
observed. 


Figure  4.7:  5-PSI.  Differential  Transducer  Calibration  Results. 

4.2.3  TBL  Measurement  Arrau  and  Data  Acquisition  System 
The  measurement  array  consisted  of  three  Endevco,  model  number 
8514-10,  pressure  transducers  (described  in  Chapter  3)  and  also,  three 
Enoevco  brand,  piezo-electric  accelerometers  (nominal  sensitivities  of  45 
mV/g).  The  array  was  located  approximately  7.5  meters  downstream  from 
the  entrance  to  the  test  section,  in  the  fully  developed,  turbulent  flow 
region  of  the  pipe. 

The  transducers  were  mounted  in  individual  cylindrical,  plug-type 
assemblies.  These  assemblies  were  inserted  into  accommodating  sleeves 
located  along  the  tunnel  wall  and  separated  circumferentially  by 
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120-degrees.  The  sleeves  were  equipped  with  gate  valves,  allowing  the 
removal  of  the  plugs  without  draining  the  tunnel.  Alignment  bolts  were 
also  positioned  on  the  outer  portion  of  each  sleeve  assembly  to  allow 
proper  alignment  of  each  plug  with  the  inner  tunnel  wall.  However,  due  to 
excessive  play  (*1  cm)  in  the  bolt  mechanism  proper,  flush  alignment  of 
the  plug  with  the  inner  tunnel  wall  could  only  be  achieved  through  tactile 
verificatioa  The  pressure  transducers  were  sealed  into  the  face  of  each 
plug  assembly  with  RTV  sealant  and  were  in  direct  communication  with  the 
flow.  Accelerometers,  were  used  to  measure  the  local  vibrations  of  each 
plug  subassembly  and  were  mounted  on  the  back  of  the  subassembly  of  each 
transducer  (see  Figure  4.8). 


Figure  4.8:  Cross-sectional  View  of  Plug  Assembly. 

The  10-volt  power  supply  used  in  the  calibration  of  the  transducers  was 
also  used  in  the  actual  tests.  Also,  the  same  cables  and  preamplifiers 
were  used  to  ensure  that  all  calibration  results  remained  unchanged. 
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All  voltage  signals—transducer  and  accelerometer— were  filtered  and 
amplified  on  Brookdeal  Differential  Amplifiers  before  processing  on  a 
Spectral  Dynamics  (SD)  375  real-time  spectral  analyzer.  The  half  power 
points  of  the  high  and  lowpass  filters  were  set  at  l  Hertz  and  1000  Hertz, 
respectfully.  The  SD  375's  analysis  bandwidth  was  set  by  the  maximum 
frequency  of  analysis.  For  example,  the  SD  375  analyzer  stores  and 
displays  400  bins  of  information;  this  sets  the  bandwidth  of  analysis  at, 

B-W-  =  Jmax  /  40°- 


Throughout  this  investigation  a  Hanning  (cos2)  analysis  window  was  used 
because  of  its  characteristic  side  lobe  suppression  (18  dB/octave)  and 
relatively  narrow  main  lobe  bandwidth.  Various  types  of  spectra  were 
chosen  from  a  selection  menu  on  the  SD  375  and  after  analyzing,  were 
either  sent  to  a  hard  copy  digital  plotter  or,  to  the  ARl/VAX  782  computer 
via  a  standard  IEEE  data  bus. 

4.3  TBL  Wall  Pressure  Spectra 

Before  the  TBL  pressure  fluctuations  were  measured  an  estimation  of 
the  vibrational  contamination,  present  in  the  difference  signals  was 
obtained.  This  estimate  was  acquired  by  measuring  the  coherent  output 
power  (COP)  spectrum  between  the  difference  signal  of  two  pressure 
sensors  and,  the  difference  signal  of  their  corresponding  accelerometer 
outputs.  As  pointed  out  in  Chapter  2,  if  the  COP  spectrum  is  significantly 
less  than  the  autospectrum  of  the  pressure  difference  signal  then  all 
vibrational ly  induced  pressure  terms  can  be  ignored.  Figure  4.9  shows  a 
typical  plot  of  the  autospectrum  of  a  pressure  difference  signal  and  the 
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corresponding  COP  spectrum.  As  Figure  4.9  indicates,  all  vibrational  ly 
induced  pressure  terms  are  essentially  removed  in  the  difference  signal. 
This  implies  that  the  power  spectrum  of  the  TBL  wail  pressure  fluctuations 
is  given  simply  by  the  cross-spectrum  between  two,  pressure  uifference 
signals  (electronic  noise  source  terms  included).  For  example, 


GAtbi(f)  =  GA“B,A-C*GAe- 
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The  turbulent  boundary  layer  (TBL)  wall  pressure  spectra  were  obtained 
as  a  function  of  glycerine  temperature  by  the  method  outlined  previously  in 
Chapter  2.  Measurement  of  the  electronic  noise  spectrum  was  found  to  be 
insignificant  (70  dB  down)  so  that. 


°Atbl^>  =  GA'B.A“C. 


Difference  signals  were  obtained  from  the  Brookdeal  Differential 
Amplifiers  and  the  corresponding  cross-spectral  density  functions  were 
obtained  from  the  SD  375  spectral  analyzer.  Figure  4.10  shows  the  TBL 
pressure  spectrum  for  various  glycerine  temperatures  where  256 
consecutive  ensemble  averages  were  taken.  These  results  indicate  that  as 
the  temperature  of  the  flow  increases,  the  frequency  content  of  the  wall 
pressure  spectrum  also  increases.  However,  the  peak  levels  of  the  power 
spectra  for  the  various  glycerine  temperatures  remain  essentially  constant, 
independent  of  Reynolds  number.  An  explanation  of  this  observation  can  be 
attributed  to  changes  in  the  viscosity  of  the  fluid  that  occurred  over  the 
temperature  range  tested.  At  the  lower  tunnel  (glycerine)  temperatures, 
the  small-scale  eddies  that  produce  the  higher  frequency  content  of  the  TBL 
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pressure  spectrum  are  more  easily  damped  by  the  fluid's  higher  viscosity. 
As  the  temperature  increases  the  viscosity  decreases  thereby  allowing 
more  of  the  smaller  scale  eddies  to  exist  which  in  turn,  generate  higher 
frequency  energy.  This  premise  is  further  substantiated  by  Figures  4.11 
and  4.12.  Both  of  these  plots  show  previously  used  non-dimensional izing 
methods  for  characterizing  the  TBL  wall  pressure  spectra  in  various  test 
mediums— air,  water  and  glycerine. 

Figure  4.11  suggests  that  the  amplitude  of  the  TBL  pressure  spectrum 
scales  with  the  inner  flow  parameter  rw,  and  outer  flow  parameters,  U  and 
pipe  diameter,  D.  The  frequency, )  is  scaled  in  the  usual  Strouhai  fashion 
where  the  length  scale  is  taken  as  the  pipe  diameter  and  velocity  as  U. 
Figure  4.12,  scales  the  frequency  in  a  similar  fashion  where  S  is  D/2. 
However,  the  amplitude  of  the  pressure  fluctuations  is  scaled  entirely  on 
outer  flow  parameters  p,  U  and  D.  Although  reasonable  data  collapse  is 
not  achieved  in  either  scaling  case,  an  interesting  result  can  be  noted.  As 
the  viscosity  (u)  of  the  medium  is  reduced,  an  increase  in  the  frequency 
content  of  the  TBL  pressure  spectrum  is  observed.  That  is.  water  having 
the  smallest  kinematic  viscosity  of  the  three  mediums  shown,  also  has  the 
highest  frequency  content  of  the  TBL  pressure  spectrum.  Glycerine  on  the 
other  hand,  has  the  lowest  frequency  content  of  the  TBL  pressure  spectrum 
and  as  expected,  the  highest  viscosity.  Also,  both  figures  suggest  that  the 
scaling  of  the  amplitude  of  the  TBL  pressure  fluctuations  is  not  critical. 
That  is,  both  methods  of  amplitude  scaling  produce  similar  results. 
However,  it  is  considerably  more  understandable  to  scale  the  amplitude  of 
the  pressure  fluctuations  with  rw,  since  this  parameter  describes  the 
physics  of  the  flow  locally,  at  the  transducer  location.  In  general,  it  is 
expected  that  the  proper  way  to  scale  the  TBL  wall  pressure 
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spectra  (both  magnitude  and  frequency)  is  entirety  with  inner  wall 
variables,  taking  due  account  of  the  resolution  imposed  by  the  finite  size 
of  the  measurement  transducer. 

In  many  previous  investigations,  the  major  parameter  used  to 
characterize  the  transducer  resolution  has  been  d/S*.  Where  d,  represents 
the  hydrophone  diameter  and  $*,  the  displacement  thickness.  Schloemer 
(1966)  recognized,  as  early  as  1966,  that  a  'small'  transducer  must  be 
small  relative  to  the  viscous  scale,  or  'wail  unit'  (o/ur),  where  ut  is  the 
friction  velocity.  This  idea  was  later  expanded  upon  by  Fabula  in  the  work 
published  by  Patrick(l977)  in  1977.  Willmarth  (1975)  also  points  out  the 
importance  of  this  form  of  scaling,  that  was  later  used  extensively  by 
Schewe  (1983).  Fabula  contended  that  the  d/S**  method  of  normalization 
fails  to  describe  the  fundamental  characteristics  of  turbulent  boundary 
layer  flow.  Fabula  argued  that  the  ultimate  transducer  resolution  is  more 
correctly  attained  by  a  maximum  value  of  urd/o  (=  d4)  rather  than  d/S**.  In 
other  words,  d4  is  a  transducer  resolution  parameter  that  is  based  on  inner 
boundary  layer  similarity  considerations,  while  d/S*  is  appropriate  for  a 
transducer  resolution  parameter  based  on  the  outer  boundary  layer  (lower 
frequencies).  In  fully  developed  turbulent  pipe  flow  S*.  is  fixed  at  roughly 
1/7 tri  ihe  pipe  radius.  From  this  it  is  possible  to  see  that  even  if  d/S* 
equals  a  constant,  d4  can  vary  considerably  over  a  given  Reynolds  number 
range  which,  varies  (with  d/S*  fixed)  from  medium  to  medium.  The 
parameter,  d4  can  be  thought  of  as  a  quantity  that  relates  the  transducer 
size  to  the  height  of  the  viscous  sublayer  region  in  the  turbulent  boundary 
layer  (y4=yur/u<5).  Fabula  proposed  that  the  ultimate  in  transducer 
resolution  is  achieved  when  urd/d  <  10.  However,  Schewe  (1983)  states 
complete  resolution  of  the  essential  structure  of  the  turbulent  pressure 
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fluctuations  with  a  measured  value  of  d+  =19,  Values  of  this  parameter 
for  data  obtained  in  the  current  investigation  ranged  from  2.1  <urd/o< 
4.4.  Therefore,  these  data  provide  the  best  spatial  resolution  of  the  Tai¬ 
wan  pressure  fluctuations  measured  to  date;  they  represent  a  true  “point" 
measurement.  (It  should  be  noted  that  the  transducer's  casing  diameter 
was  used  in  the  normalization  and  the  actual  value  of  is  most  likely 
smaller  by  a  factor  of  3,  Galib  and  Zandina  (1984)). 

Inner  wall  scaling  requires  that  the  frequency  be  non-dimensional ized  by 
the  viscous  time  scale.  Therefore,  the  frequency  of  the  wall  pressure 
spectrum  is  now  scaled  as  fu/ur2.  Figure  4.13,  while  using  the  same 
amplitude  scaling  as  in  Figure  4.11  uses  the  inner  flow,  frequency  scaling. 
This  figure  clearly  shows  the  dramatic  effect  transducer  size  has  on  the 
measured  spectra.  As  d*  decreases,  more  high-frequency  (small-scale) 
information  is  detected  by  the  transducer. 

For  completeness,  Blake*s(l970)  and  Schewe’s  (1983)  data,  obtained  for 
turbulent  flat  olate  flow,  is  compared  to  the  data  obtained  in  the  Glycerine 
tunnel  in  Figure  4.14.  In  this  figure,  the  pressure  amplitude  and  frequency 
are  scaled  entirely  on  inner  flow  parameters.  Although  internal  pipe  flows 
differ  from  external  flat  plate  flows  for  a  number  of  reasons;  the  physical 
causes  of  the  pressure  fluctuations  should  be  similar.  Therefore, 
heuristically  speaking,  if  the  parameters  that  characterize  the  pressure 
fluctuations  are  explicitly  known  then  direct  comparisons  of  the  resulting 
TBL  pressure  fluctuations  for  different  flow  geometries  is  somewhat 
justified.  At  values  of  d+>20  (Blake)  spatial  filtering  of  the  TBL  pressure 
fluctuations  is  apparently  occurring  because  of  the  observed  attenuation  of 
the  high-frequency  spectral  energy.  For  values  of  d*  <  20.  the  present  data 
and  those  of  Schewe  agree  quite  well;  attenuation  of  the  high-frequency 
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content  of  the  TBL  pressure  fluctuations  is  virtually  non-existent. 
Schewe's  contention,  that  the  complete  resolution  of  the  TBL  pressure 
fluctuations  exists  when  d+<  19  is  therefore  substantiated  by  the 
agreement  observed  between  the  two  sets  of  data.  It  appears  that  no 
further  increase  in  transducer  resolution  is  expected  for  values  of  d* 
smaller  than  20.  This  confirms  that  transducer  resolution  depends 
primarily  on  values  of  uTd/d  and  not  on  d/S*.  Scaling  the  wall  pressure 
spectra  in  the  manner  used  in  Figure  4.14  allows  the  direct  comparison  of 
various  types  of  flow  geometries;  for  example,  pipe  and  fiat  plate  flows. 
It  also  suggests  that  the  wall  pressure  fluctuations  are  mainly  controlled 
by  the  physics  of  the  flow  in  the  inner  wall  region. 

The  r.m.s.  value  of  the  TBL  pressure  fluctuations  is  obtained  by 
integrating  the  wall  pressure  spectrum  over  all  frequencies  and  represents 
a  measure  of  the  intensity  of  the  process.  Figure  4.15  illustrates  the 
dependence  of  the  r.m.s.  pressure  on  4'.  Here,  * :.r.  r.m.s.  pressure  is 
non-dimensional  ized  by  the  dynamic  head,  q<». 

Corcos  (1963)  proposed  a  correction  to  account  for  the  finite  size 
measurement  transducer  (see  Appendix  4).  This  correction  is  a  function  of 
the  non-dimensional  frequency,  u>R/Uc  and  is  applied  directly  to  the 
spectral  magnitude  of  the  TBL  pressure  fluctuations.  From  measurements 
of  the  wavenumber- frequency  response  of  a  circular  pressure  sensor, 
virtually  no  correction  is  needed  for  values  of  o)R/Uc<i.O.  From  this 
result  and  those  from  Appendix  4,  no  correction  should  be  needed  for  any  of 
the  measurements  reported  here  for  the  Glycerine  tunnel  where. 
u>max R/uc=0-6  ^max  W3S  taken  as  the  50  ^B  down  point  in  TBL  wall 
pressure  spectrum).  In  Table  4.1,  both  corrected,  and  uncorrected  values 
of  the  r.m.s.  pressure  for  Schewe's  data  is  presented  and  compared  to  those 


Figure  4. 1 5:  Dependence  of  the  Normalized  r.m.s.  Value  of  the  Wall 
Pressure  Fluctuations  upon  the  Normalized  Transducer 
Diameter. 


Investigator 

d* 

Measured 

Corrected 

Schewe 

19 

0.0098 

0.0105 

Current 

2.1 

0.0106 

0.0106 

Current 

4.4 

0.010? 

0.010? 

Table  4.1:  Comparison  of  the  Measured  and  Corrected 
Normalized  Wall  Pressure  r.m.s.  Levels. 


measured  in  the  Glycerine  tunnel.  As  indicated,  the  r.m.s.  value  of  the  wall 
pressure  fluctuations  is  well  resolved  and  equal  to  0.0107  q*,. 

4.4  Correlation  of  the  TBL  in. the  Cross-Difference  Spectrum 

The  hypothesis  that  the  TBL  pressure  fluctuations  were  uncorrelated 
between  pressure  sensors  is  a  major  assumption  of  the  TBL  measurement 
technique  described  in  Chapter  2.  This  model  assumption  may  not  be 
entirely  valid  at  very  low  frequencies  and,  if  not,  may  actually  attenuate 
the  low-frequency  part  of  the  TBL  wall  pressure  spectrum.  A  relatively 
simple  analysis  can  be  performed  to  show  that  if  the  above  assumption  is 
fully  valid  then  the  value  of  the  coherence  function  between  any  pair  of 
(pressure  sensors)  difference  signal  outputs  should  equal  0.25  over  the 
frequency  range  of  validity.  Moreover,  this  analysis  also  shows  that  the 
vibration-induced  and  acoustical  pressure  components  are  removed  in  the 
subtraction  of  any  two  pressure  sensor  signals. 

To  start,  assume  that  the  instantaneous  pressure  sensor  outputs  are 
represented  as  a(t),  b(t)  and  c(t).  Also,  assume  that  all  acoustic  and 
vibrationally  induced  pressure  signals  are  entirely  cancelled  when  any  two 
pressure  sensor  outputs  are  subtracted.  Then  using  the  same  notation  as 
in  Chapter  2.  the  Fourier  transformed  quantities  between  pressure 
difference  signals,  a(t)  -  b(t)  and  a(t)  -  c(t)  can  be  represented  as. 

A(|)-B(|)-At6|(|)-BtD,(|) 

and 


A(|)-C(|).Atol(|)-Ctbl(|). 


ElA*B)  =  E[B*C] 


0. 


Equation  4.6  will  then  reduce  to. 


ElA*A]E(A*A) 


(E(A  Al  ♦  E[B*B])(E(A  A]  ♦  ElC  C]) 


Also,  because  the  flow  is  axisymmetrical,  the  autospectra  of  the  turbulent 
boundary  layer  pressure  fluctuations  will  be  the  same  at  locations  a.  b  and 
c.  This  implies  that  (recall  that  the  “tbl“  subscript  is  suppressed). 


ElA*A]  =  E(B*Bl  =  E[C*Cl. 


Equation  4.7  then  reduces  to. 


E[A*A]ElA*Al 


(ElA  A1*EIA  AlXElA  Al*E(A  Al) 


From  this  analysis  it  can  be  stated  that  if 


1)  the  difference  signal  effectively  removes  all  acoustic  and 

vibrationally  induced  pressure  signals,  ^nd 

2)  the  pressure  sensor  separation  distance  is  large  enough  to 
permit  the  TBl  pressure  fluctuations  to  remain  uncorrelated 
between  pressure  sensors  then. 
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the  value  of  the  coherence  function  between  any  pair  of  difference  signals 
should  equal  0.25.  A  simple  measurement  will  therefore  either  prove  or, 
disprove  the  validity  of  the  TBl  measurement  assumptions. 

A  typical  measurement  of  the  coherence  function  between  two  pairs  of 
difference  sig  ,  is  shown  in  Figure  4.16  where  256  spectral  averages 
were  takea  At  very  low  frequencies  (<  10  Hert2)  the  coherence  function  is 
essentially  0.25  implying  that  all  vibrational  and  acoustic  components  of 
the  pressure  field  are  cancelled  thereby,  leaving  only  the  TBL  pressure 
fluctuations.  As  the  frequency  of  the  coherence  function  increases  above 
20  Hertz,  deviations  from  0.25  begin  to  appear.  These  deviations  are  most 
likely  due  to  the  increasing  loss  of  the  coherent  structure  in  the 
vibrational ly  induced  pressure  field  and.  not  in  the  propagating  acoustic 
field  because  the  (duct)  cutoff  frequency  for  axially  plane  wave  propagation 
is  approximately  4100  Hertz.  This  means  that  all  acoustic  noise  should 
remain  perfectly  correlated  between  pressure  sensors  for  frequencies  as 
high  as  4100  Hertz.  The  effect  of  the  vibration  contamination  is  most 
prevalent  at  lower  tunnel  temperatures  and.  for  frequencies  above  200 
Hertz,  as  shown  in  Figure  4.16.  These  vibrational  effects  appear  as 
spectral  peaks  and  valleys  in  the  TBL  pressure  spectra  (see  Figure  4.10),  but 
only  represent  a  small  fraction  of  the  total  spectral  energy. 

Another  way  of  determining  how  well  the  TBL  pressure  spectrum  was 
measured  is  by  comparing  the  power  spectrum  of  a  difference  signal  to 
cross- difference  spectrum.  From  previous  relationships  it  is  easily  shown 
hat, 


gA-bW  =  2  lim  (  £(Att)t -  B^| )*,(Att)| -Bttj) )l  /  T  =  2  6^(0  . 
T~*  00 


and  from  Chapter  2, 


6A-B,A-C$  =  Gtbl(f}- 


Therefore,  the  power  spectrum  of  a  difference  signal  should  De  3  dB  above 
the  cross-difference  spectrum  (Gj0|_).  Typically,  Figure  4.17  shows  these 
two  power  spectra.  As  expected,  the  power  spectrum  of  the  difference 
signal  A(J)-B(f).  is  nominally  3  dB  above  that  of  the  TBL  spectrum  over  the 
entire  frequency  range  showa 

The  above  analyses  and  their  results  are  significant  in  that,  vibration 
contamination  usually  causes  considerable  errors  in  the  measurement  of 
the  TBL  wall  pressure  spectrum  at  very  low  frequencies  and  has  generally 
made  accurate  conventional  measurements  (using  single  pressure 
transducers)  of  these  types  impossible  (Willmarth  (1975)).  However,  TBL 
wall  pressure  measurements  that  involve  transducer  difference  signals  can 
not  only  remove  significant  amounts  of  acoustic  and  vibratory  noise  from 
the  TBL  wall  pressure  spectrum  (Section  4.5),  but  also,  provide  a  direct 
way  to  check  the  integrity  of  the  data  acquired  (especially  at  very  low 
frequencies).  Data  integrity  of  the  TBL  wall  pressure  spectrum  can  now  be 
established  with  a  measurement  of  the  coherence  function  between 
difference  signal  pairs.  This  measurement  will  equal  0.25  for  the  precise 
measurement  of  the  TBL  wall  pressure  spectrum.  It  should  be  noted,  that 
this  difference  signal  approach  to  measure  the  TBL  wall  pressure  spectrum 
Is  not  limited  to  axisymmetric  flows;  this  type  of  measurement  technique 
will  work  as  long  as  the  transducer  array  is  positioned  normal  to  the  mean 
flow  field.  These  measurement  techniques  are  the  preferred  way  to 
measure  the  TBL  wall  pressure  spectrum. 


At  the  onset  of  this  experimental  investigation  it  was  hoped  that  the 
noise  associated  with  the  tunnel  operation  would  be  many  orders  of 
magnitude  below  that  created  by  the  TBL  wall  pressure  fluctuations. 
Because  measurements  of  the  low-wavenumber  content  of  the  TBL  have  yet 
to  be  performed  with  any  degree  certainty,  it  was  an  early  goal  of  this 
investigation  to  qualify  the  glycerine  tunnel  for  such  measurements. 
Unfortunately,  instead  of  qualifying  the  Glycerine  tunnel  as  a  candidate  for 
future  low-wavenumber  research,  measurements  of  the  acoustic  and 
vibrational  noise  spectra  have  indicated  that  these  two  sources  dominate 
the  output  of  all  flush-mounted  pressure  sensors.  Figures  4.18,  and  4.19 
show  three  different  power  spectra  measured  in  the  glycerine  tunnel,  at 
temperatures  of  35°  and  46.1°  C.  Both  plots  indicate  typical  results  of  the 
TBL  pressure  spectrum  (GttJ|),  the  acoustic  or,  cross-spectrum  (GAB)  and, 
the  power  spectrum  measured  by  a  single  transducer  (6/^.  These 
results— at  both  temperature  extremes— clearly  indicate  the  power  of  the 
cross-difference  TBL  measurement  technique  and  the  high-level  noise 
present  In  the  facility.  Over  most  of  its  frequency  range,  the  TBL  pressure 
spectrum  is  10  to  15  dB  below  that  of  the  acoustic  (and  single  transducer 
spectrum)  noise  spectrum,  and  at  46.I°C,  as  much  as  30  dB  below. 

Figures  4.20,  4.21,  4.22,  and  4.23  depict  the  corresponding  phase  angles 
of  the  TBL  pressure  spectrum  and  the  acoustic  noise  spectrum  at  glycerine 
temperatures  of  35  and  46.1°  C,  respectfully.  All  phase  plots  of  the  TBL 
pressure  fluctuations  show  essentially  zero  phase  shift  for  frequencies  up 
to  150  Hertz.  At  higher  frequencies  (>150  Hertz)  and  especially  for 
glycerine  temperatures  of  35°C.  large  deviations  from  zero  occur  in  the 
phase  of  the  TBL  spectrum. 
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These  deviations  at  this  lower  tunnel  temperature  account  for  the  more 
ragged  appearance  of  the  TBL  pressure  spectrum  (Figure  4.10)  at 
frequencies  above  150  Hertz.  (These  deviations  were  earlier  postulated  as 
being  related  to  the  loss  of  the  coherent  structure  in  the  vibrational ly 
induced  pressure  field  between  transducers  for  frequencies  above  ISO 
Hertz.)  The  plots  illustrating  the  measured  phase  difference  between 
pressure  sensors  (Figures  4.21  and  4.23)  show  that  the  relative  phase 
difference  between  pressure  sensors  is  insignificant.  This  further 
substantiates  that  all  acoustic  noise  in  the  facility  propagates  as  plane 
waves,  constant  in  phase  and  amplitude  in  any  plane  perpendicular  to  the 
flow  of  the  test  section  for  frequencies  below  1000  Hertz. 

Figure  4.24  shows  the  COP  spectrum  between  the  output  of  single 
transducer  and  its  corresponding  acceleratioa  Chapter  2  described  how 
this  particular  spectrum  represents  the  vibrational  contamination  present 
in  the  autospectrum  of  a  single  transducer.  Although  this  measurement 
indicates  a  cause  of  possible  contamination  it  does  not  indicate  the  source. 
That  is,  the  COP  spectrum  does  not  show  whether  the  noise  source  is 
entirely  structure  borne  or,  fluid/structure  coupled.  However,  Figure  4.24 
does  show  that  vibrational  effects  dominate  the  high-frequency  part  of  the 
autospectrum  of  a  single  transducer  and  that  acoustical,  fluid-borne  noise 
contributes  to  most  of  G^j)  below  150  Hertz. 

Although  not  investigated  in  detail,  it  is  speculated  that  a  major  source 
of  noise— both  acoustic  and  vibrational— is  due  to  the  tunnel's, 
pump/motor  coupling  and  the  pump.  Excessive  floor  vibration  can  be  felt  at 
the  initial  start-up  of  the  facility  and  during  tunnel  operation. 
Measurements  of  the  motor  and  pump  shaft  misalignment  were  made  and 
found  to  be  quite  excessive.  An  improved  flexible  coupling  for  the  facility 
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has  been  designed  and  subsequent  measurements  of  the  COP  spectrum  and 
transducer  cross-spectrum  should  be  repeated  after  the  new  coupling  is 
installed.  This  relatively  easy  implementation  should  reduce  tunnel 
vibrational  (and  acoustic)  noise  levels  significantly.  Other  quieting 
measures  are  also  under  investigation  at  this  writing. 

4.6  Conclusions  and  Areas  of  Future  Research 

The  techniques  employed  in  measuring  and  validating  the  TBL  wall 
pressure  spectra  were  of  primary  importance  in  this  study.  Without  these 
novel  signal-processing  techniques,  the  ability  to  measure  the  wall 
pressure  spectrum  accurately,  with  conventional  single  point 
measurements,  would  have  proved  futile  due  to  the  excessive  facility 
background  noise. 

Wall  pressure  spectra  obtained  for  both  pipe  and  flat  plate  flows  have 
been  compared  with  remarkable  agreement.  From  these  comparisons  and 
selected  scaling,  the  characteristic  features  of  the  wall  pressure  spectrum 
have  been  identified.  That  is,  the  inner  flow  parameters  rw.  and  Ur  along 
with  the  viscosity  (u)  of  the  fluid  can  be  used  to  consistently  characterize 
the  TBL  wall  pressure  fluctuations  and  their  frequency  content  for  botn 
internal  and  external  flows. 

Due  to  the  large  viscous  sublayer  present  in  the  glycerine  turbulent 
boundary  layer,  the  value  of  the  normalized  transducer  diameter  (d*)  ranged 
from  2.1  to  4.4,  depending  on  the  fluid  bulk  temperature.  These  are  the 
smallest  values  of  d*  ever  reported.  From  these  results,  a  critical 
uppermost  limit  on  the  value  of  d*  has  been  validated;  namely,  df<  20  as 
speculated  by  previous  investigators.  This  value  of  Q*  represents  an  upper 
limit  for  which  proper  resolution  of  the  high-frequency  content  of  the  TBL 
pressure  fluctuations  can  be  achieved. 


A  simple  measurement  of  the  coherence  function  between  difference 
signals  has  validated  all  the  TBL  pressure  fluctuation  measurement 
assumptions  previously  outlined  in  Chapter  2.  From  these  results, 
measurement  of  the  very  low-frequency  content  (<10  Hertz)  of  the  TBL 
wall  pressure  spectrum  was  obtained.  The  coherence  function 
measurement,  also  provides  a  simple  means  to  test  the  integrity  of  the  TBL 
wall  pressure  spectrum  (i.e.  removal  of  vibration  and  acoustic  noise), 
especially  at  very  low  frequencies. 

An  area  of  possible  future  research  for  the  glycerine  tunnel  may  include 
velocity  and  pressure  correlation  measurements.  Because  the  wall 
pressure  fluctuations  are  related  to  the  fluctuating  velocity  field  and. 
because  the  viscous  sublayer  is  quite  thick,  detailed  measurements  of  these 
types  are  possible.  They  would  provide  additional  information  on  coherent 
structures  including  bursting  frequencies  and  the  scaling  of  these 
frequencies.  A  major  issue  in  the  study  of  coherent  structures  in  turbulent 
boundary  layers  has  been  the  scaling  of  the  bursting  frequency.  This 
controversy  centers  about  transducer  spatial  resolution.  It  has  been  shown 
that  (conservative)  d*  values  as  small  as  2.1  can  be  obtained  in  the 
glycerine  tunnel,  which  would  therefore  make  the  glycerine  tunnel  a  likely 
candidate  to  resolve  this  controversy. 


APPENDIX  1 

SPECTRAL  DENSITY  VIA  THE  FINITE  FOURIER  TRANSFORM 


The  two  sided  spectral  density  functions  $xx(j),  and  SXy(f)  are  by 
definition,  the  Fourier  transform  of  their  auxo  and  cross-correlation 
functions  Rxx(r).  and  RXy(r).  However,  they  can  also  be  represented  as  the 
expected  value  of  their  Fourier  transformed  time  functions  x(t)  and  y(t). 
Before  proceeding  with  this  proof  a  few  fundamental  definitions  will  be 
given 

A  random  process  can  be  described  as  an  ensemble  of  a  real-valued  time 
function  and  is  denoted  as  {xk(t)}  where,  -«o  <  t  <  «>.  The  function,  xk(t)  is 
defined  as  the  kth  event  of  an  entire  ensemble  in  which  N  events  occur  or 
are  observed.  In  practice,  N  is  a  statistically  large  finite  number  and  is 
usually  large  enough  to  approximate  infinity.  The  mean  values  of  the 
stationary  random  processes  (xk(t)}  and  {yk(t)},  are  described  as: 


00 


jix  =  E(xk(t)]  =  j  x(t)  p(x)  dx 


-00 


00 


Jiy  =  E[yk(t)l  =  Jy(t)  p(y)  dy 


-00 


where  p(x)  and  p(y)  are  the  probability  density  functions  of  x(t)  and  y(t), 
respectfully.  For  arbitrary  fixed  t  and  z,  let: 


Ryy(r)  =  E[yk(t)  yk(t  ♦  t)l 

Rxu(t)  =  ECxk(t)  yk(t  ♦  t)i 


Now,  consider  the  kth  sample  records  xk(t)  and  yk(t)  of  stationary  random 
processes  (xk(t)}  and  (yk(t)}.  Let, 

SxyCf.TJc)  =  Xic^Ci.T)  Yk(f.T)/T  (A.l) 

where,  0  <  t<  T  and, 


Xk(j ,T)  =  J  xk(t)  e'J2lTfl  (A-2a) 


Yk(| J)  =  j  yk(t)  e' 
0 


(A.2.0) 


The  above  integral  equations  represent  the  finite  Fourier  transform  of  the 
records  xk(t)  and  yk(t),  respectfully.  We  now  wish  to  show  that. 


Sxy(f)  =  lim  E[  Sxy(I,T,k)  ]  /  T. 

J  ->oo 


Substituting  Equations  A.2  into  Equation  A.l  and  changing  variables  to 
avoid  confusion  results  in, 
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T  T 

Sxy(J.T,k)  -  '/t{  J  xk(0  yk(0  e')2") ‘  ^  dj  OC-  (Al.3) 

0  0 

Now,  let  r  =  OS  and  dr  =  dC  so  the  range  of  integration  changes  from  (CO 
to  (£,r).  Equation  Al.3  now  becomes, 

0  T 

Sxy(J,T.k)  =  }  [’  /t  j  xk(0  yk(Or)  dS  ]  e‘J2lTlr  dr 
-T  -r 

T  T-r 

*  j  [ 1  /t  \  xk(0  yk(Or)  d£  ]  2$^  dr.  (A1.4) 
0  0 

Talcing  the  expected  value  of  both  sides  of  Equation  A1.4  and  recalling  that 
the  definition  of  the  cross-correlation  function  is, 

Rxy(r)  =  E(xk(t)  yk(t  ♦  r)l 


Equation  A1.4  becomes, 

0  T 

E[SXy  (f.T.k)l  =  jf/.j.  xy(t)  ]  <f  l2ltf*  dr 

-T  -r 


T  T-r 

*1  ['/T  }Rxy^)dC  1  e"i27Tlrdr 
0  0 


*»*4*t^*  <  .4^*1  *.-•  «.■«  *_-«  *_r*  M  -"4. 


»  \ ( 1  -  Itl  /  T)  Rxy(r)  e‘i2TrJr  dr.  (ai.S) 

-T 

Taking  the  limit  as  T  tends  to  infinity  of  Equation  A1.5  results  in  the 
following  expression, 

00 

lim  E[5Xy(J.T,k)]  -  \  Rxy(r)  e'i27Tir  dr.  (AI.6) 

T"*00  -CO 

But  from  the  definition  of  Sxy(f,T,k),  Equation  AI.6  becomes. 

00 

s^d)  =  lim  E(Xk«(i.T)  Yk(f,T)/T|  *  } R„y(r)  e‘i2TTlr  dr.  (AI.7) 


The  above  expression  Is  the  desired  result  for  the  cross-spectral  density 
function  The  other  auto  spectral  density  functions  of  x(t)  and  y(t)  follow 
directly  and  are  given  by, 

00 

Sxx(f)  =  lim  E[Xk"(f,T)  Xk(f,T)/Tl »  \ Rxx(r)  e")2TTfT  dr  (A1.8) 


Syyd)  =  lim  ElYk-(f.T)  Yk(J.T)/T) »  jRyy(r)  e~)2nlz  dr.  (AI.9) 


-CO 


to  and  cross-spectra?  density 


l<x(r)e’j2TTJt  dr  (AUO) 


,y(r)  e~i2TTlr  dr  (Al.ll) 
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6XU(|)  =  2  lim  ElXk"(j.T)  Yk(|,T)/T)  -  2  {fUtr)  e'J2TTfr  dr.  (A1.12) 


The  above  expressions  show  the  relationship  between  the  cross  and  auto 
spectral  density  functions  and  the  Fourier  transforms  of  their  time  records 
xk(t)  and  yk(t).  These  expressions  are  exact  and  represent  the  true  spectral 
functions  for  the  limiting  case  of  T. 


APPENDIX  2 

SOUND  PROPAGATION  IN  CIRCULAR  TUBES 


The  general  wave  equation  in  cylindrical  coordinates  is  given  as. 

a2*  i  a*  i  a2*  a2*  i  a2* 

- ♦ - + - ♦  -  = - (A2.1) 

3r2  r  3r  r2  8<p2  922  c2  at2 

In  the  above  expression  r,<p .  and  z  are  the  usual  spatial  coordinates,  t  is  the 
time  variable  and  c.  the  sound  speed  of  the  medium.  Solutions  to  Equation 
A2.I  will  be  of  the  form. 

*  =  R(r)2(z)W)T(t).  (A2.2) 

Substituting  equation  A2.2  into  Equation  A2.I  and  using  separation  of 
variables;  the  wave  equation  for  the  radial  dependency  is  given  as. 


1  82R  1  3R 

R  9r2  rR  dr 


kr2  -  m2/r2  =  0. 


(A2.3) 


In  the  above  expression  kr.  and  m  are  separation  constants.  The  standing 
wave  solution  to  Equation  A2.3  is  given  in  terms  of  Bessel  functions  and 
Neumann  functions  of  order  m.  such  that; 

R(r)  =  A  Jm(kr  r)  *  B  Nm(  kr  r).  (A2.4) 

where  A  and  B  are  constants,  for  wave  propagation  in  a  cylinder  of  outer 
radius  a,  the  solution  must  be  defined  at  r=0,  therefore  B  must  be 


identically  equal  to  0.  For  a  tube  with  rigid  walls  the  velocity  is  2ero  at 
r=a.  Therefore,  the  roots  of  Equation  A2.4  are  computed  from  (where 
=  kra), 

3  f)  !  =  kf  Jffl  =  0.  (A2.5) 

—  'r-a 

3r 

Table  A2.I  shows  the  first  few  roots  of  Equation  A2.5. 


Table  A2.1:  Roots  of 


m 

o=0 

o=l 

0=2 

0=3 

0 

0 

3.83 

7.02 

10.17 

1 

1.84 

5.33 

8.54 

11.71 

2 

3.05 

6.71 

9.97 

13.71 

3 

4.20 

8.02 

11.35 

14.59 

For  frequencies  below  the  lowest  non-zero  mode  (0,1)  of  propagation,  only 
axial  plane  waves  can  propagate  in  the  tube.  Therefore,  an  expression  for 
this  lowest  cutoff  mode  can  be  written  as, 


2f)o  =  a  /  c  =  1.84 


where  c,  is  the  sound  speed  of  the  fluid  and  is  the  radian  frequency.  In 
terms  of  the  lowest  cutoff  frequency,  the  above  expression  can  be  written 


The  above  expression  represents  the  uppermost  frequency  for  which  plane 
axial  wave  propagation  exists  in  a  rigid-walled  tube.  Higher  order  modes 
represent  non-axial  wave  propagation  and  are  excited  at  the  tube 
termination  or,  from  disturbances  in  the  tube. 


APPENDIX  3 

RANDOM  AND  BIAS  ERRORS  IN  THE  CROSS-DIFFERENCE  SIGNAL  SPECTRUM 
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Essentially  two  types  of  errors  occur  in  the  analysis  of  random  data. 
These  types  of  errors  are  normally  termed  the  “random*  and  “bias"  errors 
of  measurement. 

Bias  errors  (%)  are  independent  of  the  data  under  analysis  and  result 
from  the  number  of  ways  that  data  is  analyzed  and  processed.  For  instance, 
data  tapering,  a  process  that  involves  various  types  of  data  weighting,  has 
a  number  of  errors  associated  with  it.  A  major  type  of  one  of  these  errors 
(in  the  measured  data)  is  due  to  the  bandwidth  associated  with  the  analysis 
window.  That  is.  the  amplitude  of  a  harmonic  estimate  is  biased  due  to  the 
broad-band  noise  included  in  the  bandwidth  of  the  analysis  window's  major 
lobe.  Side  lobes— associated  with  most  types  of  analysis  windows— also 
allow  spectral  leakage  to  "spill  over*  into  the  major  lobe  and  can  account 
for  significant  errors  in  the  harmonic  estimate  if  large  amounts  of  higher 
frequency  energy  is  present.  These  types  of  errors  result  in  an  constant 
increase  (bias)  in  the  true  spectrum  that  would  otherwise  be  observed  if 
the  analysis  window  had  infinite  resolution  (i.e.  zero  bandwidth  main  lobe 
with  no  side  lobes). 

Random  errors  (€r)  occur  because  averaging  operations  must  occur  over 
a  finite  extent  or,  for  a  limited  sample  size.  The  random  scatter  in  the 
measured  values  that  results  from  the  finite  sampling  time  are  then 
referred  to  as  'random*  error. 

Bendat  and  Piersoi  (1980)  have  derived  various  formulae  that  describe 
the  bias  and  random  errors  inherent  in  the  auto  and  cross-spectral  density 
estimates  and  for  the  coherence  function.  These  errors  are  tabulated 
below  in  Table  A3.I. 


03 


Table  A3.1:  Estimation  Errors  in  the  fleasured  Spectral  Density  Functions. 


Estimated 

Function 

Bias  Error, 

eb 

Random  Error, 

€r 

Autospectrum, 

6xxW 

-1/3  (Be/Br)2 

1 

Cross-spectrum, 

6xy(i) 

-1/3  (Be/Br)2 

'/l?(xy(l)l/nd 

Coherence, 

V() 

Undefined 

/2:,'^y'/iv^d 

In  the  above  table  Bg  and  Br  refer  to  the  spectral  resolution  and  half 
power  point  bandwidths,  respectfully  and  r\j  is  the  number  of  spectral 
averages  taken. 

For  data  acquired  in  the  glycerine  tunnel  random  errors  in  the  estimate 
of  the  cross-difference  signal  spectrum  (6^|)  are  approximately  on  the 
order  of  12$  or,  U  dB.  The  corresponding  bias  error  due  to  the  2.5  Hertz 
analysis  bandwidth  (Be)  is  negligible  if  the  turbulent  boundary  layer 
pressure  spectrum  is  approximated  as  band-limited  white  noise  (i.e.  Br  in 
this  case  will  be  the  bandwidth  of  the  data).  The  random  error  associated 
with  the  coherence  function  of  the  (pressure)  transducer  difference  signals 
is  estimated  at  135*5  or.  1.2  dB.  The  bias  error  associated  with  a  typical 


v'Vv' 

/*.•**/• 

•A'.V 


m 


TTv 

v  s  * 
/  V  v 
V  V’ 


autospectrum  measured  by  a  single  transducer  is  difficult  to  estimate  due 
to  the  sharp  peaks  tlv  occur  in  spectrum  (i.e.  hard  to  define  Br).  However, 
the  random  error  a*  ociated  with  the  autospectrum's  256  spectral  averages 
is  on  the  order  of  6*  or,  0.5  dB.  Random  errors  associated  with  the 
cross-spectral  (acoustic)  measurements  of  two  pressure  sensors  are  not 
known,  because  a  corresponding  coherence  function  measurement  was  not 
taken.  However,  these  random  errors  are  not  expected  to  be  any  greater 
than  those  associated  with  other  spectral  measurements  and  are  at  most 
between  10  to  12*. 


APPENDIX  4 

CORCOS1  CORRECTIONS  TO  THE  ESTIMATED  TBL  WALL  PRESSURE  SPECTRUM 


Cor  cos  (1963)  proposed  various  correction  factors,  due  to  the 
attenuation  of  the  high-frequency  content  of  the  turbulent  boundary  layer 
pressure  spectrum  (TBL)  caused  by  its  measurement  with  a  finite-sized 
sensor.  These  correction  are  based  on  the  similarity  parameter  a>a/Uc  in 
which  a  is  the  transducer  radius  and  Uc,  the  convection  velocity  of  the 
pressure  producing  eddies.  The  following  summarizes  Corcos'  results  for  a 
circular  pressure  sensor. 

Table  A4.I:  Attenuation  Factors  for  a  Round  Transducer. 


o>a/Uc 

W* 

10log(*  m/*) 

0.05 

0.9651 

-0.15 

0.10 

0.9313 

-0.31 

0.15 

0.8996 

-0.46 

0.20 

0.8698 

-0.61 

0.30 

0.8123 

-0.90 

0.40 

0.7585 

-1.20 

0.50 

0.7069 

-1.50 

0.60 

0.6573 

-1.82 

0.70 

0.6094 

-2.15 

0.80 

0.5632 

-2.49 

0.90 

0.5166 

-2.85 

1.00 

0.4758 

-3.23 

In  the  above  table  4>m  is  the  measured  TBL  wall  pressure  spectrum  with  a 
transducer  of  radius  a  and,  $  is  the  true  spectrum,  measured  by  a  point 
transducer. 


The  frequency  range  of  the  TBL  wall  pressure  spectrum  measured  in  the 
glycerine  tunnel  was  from  l  Hz,  to  800  Hz  (50  dB  down  point).  This  range, 
corresponded  to  normalized  transducer  radii  (<oa/Uc)  of  7.9*I0“4  to.  0.63. 
From  these  values  and  the  corresponding  corrections  given  in  Table  A4.1, 
corrections  to  the  TBI  wall  pressure  spectrum  are  for  the  most  part  less 
than  idB  and  therefore  insignificant,  considering  the  accuracy  of 
measurement.  In  closing,  if  the  true  parameter  describing  the  high 
frequency  resolution  of  the  TBL  wall  pressure  spectrum  is  <f(=  urd/u),  as 
Chapter  4  suggests,  then,  the  above  corrections  based  on  the  similarity 
parameter  <oa/Uc  should  not  apply. 
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